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Abstract. The structural and nucleosynthetic evolution 
of 3, 4, 5, 6 and 7Mq stars with two metallicities {Z ~ 
0.005 and 0.02) has been computed in detail, from the 
early pre-main sequence phase up to the thermally pulsing 
(TP) AGB phase or the onset of off-center carbon burn- 
ing. Typically 10 to 20 thermal pulses have been followed 
for each TP-AGB object. This homogeneous and quite 
large set of models allows us to present an overview of the 
thermal pulse properties as well as of the nucleosynthesis 
accompanying the TP-AGB phase of intermediate-mass 
stars. 

More specifically, after a brief description of the pre- 
vious evolutionary stages, predictions are given for the 
isotopic ratios involving C, N, O, Ne, Mg, Al and Si. 
Also the surface abundances of ^Li, ^^F and ^'^Na are re- 
ported. As the asymptotic phase of the thermal pulses 
has been reached for each star, we also indicate how these 
abundances will probably evolve until the stars completely 
loose their envelope, by including the evolution of the nu- 
cleosynthesis itself. 

This article, in its paper form, has been shortened at 
a level of roughly 60 % as required by directives coming 
from the A&A editors. The complete article (50 pages 
containing 37 figures) is only available in electronic form. 

Key v^rords: stars: AGB - stars: structure - stars: nucle- 
osynthesis - stars: abundances 



1. Introduction 

At the end of their life, stars of low- and intermediate- 
mass evolve along the Asymptotic Giant Branch, where 
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they experience recurrent thermal instabilities and sub- 
stantial mass loss. During this phase, the stars undergo a 
very rich and unique nucleosynthesis. Moreover, recurrent 
dredge-up events enrich the stellar surface with the freshly 
synthesized nuclides which are then ejected into the inter- 
stellar medium through the strong winds. The thermally 
pulsing AGB stars thus play a crucial role in the chemical 
evolution of galaxies. 

On the other hand, AGB stars are important contrib- 
utors to the integrated luminosity of stellar systems in a 
wide range of ages. A coverage of this phase is needed in 
particular for the treatment of the luminosity functions 
and color magnitude diagrams of star clusters. Last but 
not least, the corresponding stellar models are basic tools 
to study carbon stars, OH/IR stars, planetary nebulae and 
their central objects (white dwarfs). 

In the last 15 years, extensive observational data has 
become available that changed and highly constrained the 
theoretical view of the AGB evolution. Most of the pro- 
gresses came from Magellanic Cloud studies, which re- 
vealed for example the absence of very luminous carbon 
stars (which were previously predicted to originate from 
stars with initial masses between 5 and 8 Mq), and the un- 
expected existence of relatively faint and lower mass car- 
bon stars (Blanco et al. 1980, Mould & Aaronson 1982, 
1986; Wood et al. 1983; Reid & Mould 1984; Aaronson 
& Mould 1985; Smith & Lambert 1989). This so-called 
"carbon star mystery" (Iben 1981) has led to the concept 
of hot-bottom convective envelope burning (Iben 1975; 
Sackmann et al. 1974; Scalo et al. 1975; Sackmann & 
Boothroyd 1992). This mechanism, which reconverts C- 
rich to 0-rich envelopes and produces ^Li, was confirmed 
by the existence of lithium enriched AGB stars (Smith & 
Lambert 1989, 1990b; Plez et al. 1993; Smith et al. 1995). 
On the other hand, the discovery, inside meteorites, of in- 
terstellar grains which have been formed in the wind of 
cool carbon stars (see e.g. Zinner et al. 1991a), comple- 
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mcntcd the stellar abundance analysis, and brought cru- 
cial informations on AGB star nucleosynthesis. In spite of 
important improvements, many important issues of AGB 
evolution still remain to be addressed, observationally as 
well as theoretically. 

This paper is the first of a series aimed at exam- 
ining all the physical processes that may induce abun- 
dance anomalies in the surface of evolved AGB stars. In 
the present work, we specifically investigate the case of 
intermediate-mass stars. Stars with initial masses between 
3 and 7 Mq and with metallicities Z of 0.02 and 0.005 have 
been evolved, starting from the pre-main sequence until 
they had gone through a number of helium-shell flashes 
along the TP- AGB. The 7Mq (6 and 7Mq) stars with 
Z = 0.02 (0.005) ignite off-center carbon burning and will 
not be discussed in detail. For each other stellar mass, 
full evolutionary models were computed up to the fourth 
"full amplitude" thermal pulse in the so-called "asymp- 
totic regime". Synthetic evolution was then used to fol- 
low global surface properties and nucleosynthesis up to 
the planetary nebula ejection. This self- consistent and ex- 
tended set of models is particularly relevant for population 
synthesis and chemical evolution models. Very complete 
figures and tables containing informations about the cen- 
ter and surface evolution of the structure and chemical 
composition for all our models are available. 

In the present computations, which include the lat- 
est input physics, the evolution of the chemical compo- 
sition is followed by only taking into account nuclear re- 
actions and mixing inside convective zones (we assume 
instantaneous mixing, except when hot-bottom burning 
occurs, which requires a time-dependent convective diffu- 
sion algorithm). We neglect all other mixing mechanism 
of non-convective origin which may contribute to trans- 
port chemical elements in radiative regions, like rotation- 
induced mixing or diffusion. In the stellar mass range we 
consider, this approximation is supported by the observed 
surface abundances at different evolutionary stages. This 
is however certainly not the case for lower mass stars, for 
which sonic! "non-standard" mixing processes have to be 
invoked to explain various chemical anomalies. 

The physical ingredients of our models and numerical 
aspects of the computations are presented in Sects. 2 and 
3, respectively. The general properties of the stars and 
their chemical evolution during the phases prior to the TP- 
AGB one are summarized in Sects. 4 and 5, respectively. 

In Sect. 6, we briefly present the structural features of 
the thermal pulses and discuss the occurrence of the third 
dredge-up in our models. We derive a core-mass luminosity 
relation for our less massive stars that do not undergo 
hot-bottom burning. We then present the extrapolation 
procedure that we use to follow both the evolution and the 
nucleosynthesis in the "asymptotic regime" up to the TP- 
AGB tip. We also give the initial mass-final mass relation 
resulting from our models, that appears to be in good 
agreement with the observed one. 



In Sect. 7, we review the various nucleosynthesis sites 
inside TP-AGB stars, namely the burning shells (He, ^'^C 
and H), the convective tongue that develops during a ther- 
mal runaway and the base of hot convective envelopes. We 
investigate in detail the nucleosynthetic processes that in- 
volve the Li, CNO, F, NeNa and MgAl elements and neu- 
trons. We then discuss which specific nuclear region, in- 
side a TP-AGB star, mainly contributes to the surface 
abundance change of each nuclide when a third dredge- up 
occurs. 

In Sect. 8, we finally present our predictions concerning 
(i) the evolution of all the surface isotopic ratios by includ- 
ing our evolutionary models and synthetic extrapolations 
up to the convective envelope removal and (ii) the chemi- 
cal yields ejected at different ages for all the species up to 
^*Si. As in the previous sections, comparisons between our 
predictions and observations are attentively performed. 

In this paper version. Sects. 2, 3, 5 and 8 have been 
somewhat reduced. 

2. Physics of the models 

2.1. Structure 

2.1.1. Equation of state 

Our equation of state is analytic. It includes 

— the ionization of H, He, C, N and O. Complete ioniza- 
tion is assumed when OIX/OVIII > 50; 

— the electrostatic corrections through the statistical 
model of Debye-Hiickel for partially ionized regions 
and through the interpolation formalism of Koester 
(1976), between the Debye-Hiickel and the Thomas- 
Fermi models, for the completely ionized matter; 

— the electronic degeneracy as soon as its contribution 
to the total pressure exceeds 0.1%; 

— the production of neutrino pairs through four different 
processes, namely the photo-neutrinos, plasma neu- 
trinos, (e+e") pair annihilations and bremsstrahlung 
neutrinos. Energy loss contributions comes from Mu- 
nakata et al. (1985), Itoh et al. (1992) and Kohyama 
et al. (1993). 

2.1.2. Opacities 

At low temperatures, i.e. below 8000 K, we use the Alexan- 
der & Ferguson (1994) opacity tables, very well suited for 
cool red giant envelopes and atmospheres. 

At temperatures above 8000 K, we use the OPAL 
opacity tables computed by Rogers & Iglesias (1992). 
The abundance distribution of the heavy elements be- 
ing scaled from the solar composition given by Anders & 
Grevesse (1989) and updated by Grevesse (1991), as for 
low-temperature opacities. 

For both sets of opacity tables, we first interpolate bi- 
linearly in T and R = log(/!)/Tg ) for each table the compo- 
sition of which is close to the model composition and then 
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linearly in X, Z (and -'^(O)) to obtain the right opacity 
coefficient corresponding to each shell. 

Finally, we used the Hubbard & Lampe (1969) pro- 
gram to generate conductive opacity tables corresponding 
to the same chemical compositions and (i?, T) grid points 
as for the OPAL tables. 

2.1.3. Atmosphere treatment 

The stellar structure is integrated from the center to a 
very low optical depth (r — 0.001) in the atmosphere. We 
constraint the atmospheric temperature profile to corre- 
spond to those coming from realistic atmosphere models 
obtained by integrating the radiative transfer equation. 
More specifically, we use atmosphere models (1) from Plez 
(1992) for Tcft < 3900 K, (2) from Eriksson (1994, private 
communication, using a physics similar to Bell et al. 1976) 
up to 5500 K and (3) computed with the Kurucz atmo- 
sphere program above 5500 K. At each time step, we lin- 
early interpolate a T(r) profile corresponding to our model 
Teff , log g and Z from these grids of atmosphere models 
and then correspondingly modify the radiative pressure 
and gradient. 

2.1.4. Convection 

The structure of the convective regions is computed using 
the classical Mixing- Length Theory (MLT), as prescribed 
by Kippenhahn et al. (1968). Our models are standard in 
the sense that (1) the Schwarzschild criterion is consid- 
ered to delimit the convective zones and (2) neither over- 
shooting nor semi-convection have been considered in the 
present computations. The ratio a of the mixing-length 
free parameter over the pressure scale height has been put 
to a value of 1.5, i.e. rather close to qiq = 1.64 with which 
we fit the solar structure. 

2.1.5. Mass Loss rates 

^Prom the Main Sequence up to the thermally pulsing 
AGB phase, we used the empirical Reimers (1975) formula 

M = -3.981O-"r;:^M0yr-i, (1) 

with r] = 0.5 up to the central He exhaustion. 

^From the beginning of the AGB phase, we adopted 
variable rj values, depending on the total mass of the star. 
More specifically, we took -q = 2.5, 3, 3.5 or 4 for our 3, 4, 5 
or 6 Mq models respectively, whatever Z. This r] increase 
with the total mass was found by Bryan et al. (1990) to 
give rather good results along the AGB. However, later 
during the thermally pulsing phase, Blocker (1995) re- 
cently showed that the Bryan et al. (1990) formula un- 
derestimates the observed mass loss rate. This leads us 
to further increase r] (see Sect. 6.3). However, as men- 
tioned by Vassiliadis & Wood (1993), whatever empirical 
formula, the mass loss rate increase used in the models 



along the thermally pulsing AGB phase still remain un- 
certain by more than a factor of two. As demonstrated 
in Sect. 8, this important uncertainty on the mass loss 
rate formula conditions the AGB phase duration and con- 
sequently, the surface composition evolution through the 
successive third dredge-up events. 

2.1.6. Rate of gravitational energy change 

We note that along the thermally pulsing AGB phase, 
the thermal pulse features slightly depend on the way one 
estimates the rate Sgrav of gravitational energy release or 
gain (see e.g. Kippenhahn & Weigert 1991). So, let us 
specify that we wrote 

__Du P_Dp 

Sgrav - + ^2 ' 

where u is the specific (i.e. by unit mass) internal energy. 

2.2. Nucleosynthesis 

2.2.1. Selected nuclides and network 

We follow the abundance evolution, through the whole 

stars, of 45 nuclides, namely the neutrons, all the 31 stable 
nuchdes up to as well as ^Be, ^B, "C, "C, ^^N, i^q, 

ISp, 22Na, 26Alg, 26Aim^ 27^^^ 30p ^nd 31S. 

These nuclides interact through a network containing 
172 nuclear (neutron, proton and "'He captures) and de- 
cay reactions. Many nuclear reaction rates are taken from 
Caughlan & Fowler (1988). Let us just mention the nuclear 
reactions for which we selected another rate prescription. 

- i2C(n,7) from Nagai et al. (1991), 

- ^2c(q:,7) from Caughlan et al. (1985), 

- "C(n,7) from Raman et al. (1990), 

- "C (n , 7) from Beer et al. (1991), 

- "C (p,7) from Wiesher et al. (1990), 

- ^''C (a ,7) from Funck & Langanke (1989), 

- "C (a , n) , (n , 7) , (n , 7) and ^^O (n , a) 
from Fowler et al. (1967), 

- i-^N (n , p) from Brehm et al. (1988), 

- i^N {a , 7) from De Oliveira (1995), 

- 1^0 (n,7) from Allen et al. (1971), 

- ^^O (n , 7) from Wagoner (1969), 

- 1^0 (p, 7) and ^'^0(p,a) from Landre et al. (1990), 

- ^^0(n,7) from Thielemann (1991), 

- ^^F (n , 7) and ^iJSfe (n , 7) from Bao & Kappeler 

(1987) , 

- i9F(p,a) from Kious (1990), 

-20Nc(n.7), 2^Ne(n,7), 23Na(n,7), 24Mg(n,7), 
^^Mg (n , 7) and ^^Mg (n , 7) from Beer et al. (1991), 

- 22^3(0;, 7) and 22Ne(a,n) from Hammer (1991), 

- 22Na (n , 7) , 23Na (Q! , p) , 2^A18 (n , 7) and 
^^Al™ (n , 7) from Woosley et al. (1978), 

- ^%lg (a , p) and ^''Al (p , a) from Champagne et al. 

(1988) , 
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- 25Mg(p,7) and 26Mg(p,7) from Illiadis (1990), 
-and 26Als(p,7) and ^^A\B{j)'^^Ar' from Vogelaar 
(1989). 

On the other hand, the nuclear reaction rates 

for i^F (n , p) , i^F (n , a) , "^^Mg (0,7), ^SMg (a , p) , 
26]\/[g (a , 7) , ^^Mg (a , p) , ^®Mg (a , n) , (a , n) as 
well as the 45 nuclear reactions involving Si, P and S are 
fitted with a procedure developed by Rayct (1993, private 
communication). Let us finally mention recent improve- 
ments not yet incorporated in our network: various (n, 7) 
reactions rates on light nuclei (Nagai ct al. 1995) and a 
capture rates on by Kappeler et al. (1994). 

Finally, the nuclear screening factors are parameter- 
ized by using the Graboskc ct al. (1973) formalism, in- 
cluding weak, intermediate and strong screening cases. 

2.2.2. Neutron abundance 

Neutrons can be produced during the thermally pulsing 
AGB phase. They can, in particular, be captured by nu- 
clides heavier than ^^S, mainly through (n, 7) reactions. 
As we do not follow the abundance of such elements, we 
created an additional nuclide, called ^"^S^, the mass frac- 
tion of which is set to 



209 



(3) 



a=34 



This sum counts in fact 253 stable nuclides up to ^"^Bi. 
In order to have a good estimation of the free neutron 
abmidancc, we follow a method described by Jorissen & 
Arnould (1989). 

2.2.3. Mixing inside convective zones 

If nuclear reactions occur inside a convective zone, most of 
them proceed in general with a longer time scale than the 
turn-over time scale associated to the motion of convective 
cells. In such a case, one usually consider the convective 
mixing as instantaneous. This allows to treat nucleosyn- 
thesis in one shot, by taking mass- weighted averages of 
the number abundances and nuclear reaction rates over 
the whole convective region. 

However, if some key nuclides nuclear ly evolve more 
rapidly than they are mixed, we have to consider the con- 
vection transport together with the nuclear reactions, so 
that nucleosynthesis equations become difi^usion equations 
of the form 



dt 



J_d_ 
pr'^ dr 



pr DcoNV-^ 
or 



nuclear terms . 



(4) 



where Yi = Xi/Ai, the ratio of the mass fraction over the 
atomic mass of each nuclide i, and Dconv = v'^tconv, 
the convective diffusion coefficient, with v and tconv be- 
ing the mean velocity and turn-over time of the convec- 
tive cells, respectively. It is very CPU time consuming to 



treat the nucleosynthesis and convective mixing together 
through such diffusion equations for each nuclide in each 
convective shell for all the time steps. Consequently, as 
usual, instantaneous mixing has been assumed by default 
for our stellar evolution computations, except in two cases 
that potentially require a time-dependent treatment of the 
convective mixing: (1) the convective tongue associated 
with thermal pulses and (2) the base of the convective en- 
velope when it is hot enough, both along the AGB phase. 

We verified that inside thermal pulses, the elements 
that are the most affected by a diffusion treatment of 
the convective mixing are the very unstable nuclides (like 
^^Na), neutrons and protons. However, this do not change 
very significantly the nucleosynthesis global results (see 
however Sect. 7.5). 

On the contrary, the mean turn-over time scale of the 
convective motions is much longer inside the convective 
envelope than inside a thermal pulse. If nuclear burning 
occiirs at the base of the convective; c^nvelope, some im- 
portant nuclear reactions can, depending on the bottom 
temperature, occur faster than the convective mixing so 
that a diffusion treatment is required. This is in particu- 
lar the case for some of the reactions involved in the ''Li 
synthesis (see Sect. 7.2). 

2.3. Initial models 

We have built initial models of 3, 4, 5, 6 and 7Mq allow- 
ing us to start our computations at the beginning of the 
pre-main sequence phase (along the Hayashi track), when 
central temperatures are below 5 10^ K. 

For our Z = 0.005 (0.02) models, the hydrogen mass 
fraction has been set to X — 0.745 (0.687) and that of 
helium to F = 0,250 (0.293), The abundances of the heav- 
ier elements have been scaled on the Anders & Grevesse 
(1989) abundance distribution. 

3. Numerical aspects 

We solve the five one-dimensional hydrodynamical equa- 
tions for the evolution of the stellar structure with the 
mass coordinate as independent variable. 

3.1. Structure 

We first convert the five non-linear stellar structure dif- 
fca-ential equations into difference equations following the 
Henyey's scheme by defining discrete shells in mass. We 
then solve the resulting linearized system of equations by 
the Newton-Raphson relaxation procedure (see e.g. Press 
et al. 1986). Typically, 600 shells are required for a model 
on the main sequence phase, while along the AGB, 1000 
to 1300 shells are needed (200 to 500 more shells being 
added during a thermal pulse). 

On the other hand, typically 200 (700) models are nec- 
essary to model the central hydrogen (helium) burning 
phase. During the AGB phase, 2500 to 4500 models are 
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required to model the time separating two successive ther- 
mal pulses, while the thermal pulse itself requires 100 to 
350 models. 

Let us stress that these rather small numbers of mass 
shells and time steps have been chosen in order to save 
computer time and present a very large set of nucleosyn- 
thesis predictions. Straniero et al. (1996) very recently 
showed that this is probably in the disadvantage of the oc- 
currence of the third dredge-up along the TP- AGB phase. 

Finally, the mass AAf = MAt lost during each time 
step At is suppressed proportionally to the shell masses 
of the whole envelope above the burning regions. 

3.2. Nucleosynthesis 

The nucleosynthesis equations are solved using the Wag- 
oner (1969) numerical technique, well suited for our pur- 
poses. However, if a diffusion treatment is required to 
simultaneously compute the time-dependent convective 
mixing and the nuclear burning, the diffusion Eqs. (4) arc 
then solved with the same Newton-Raphson algorithm. 

4. Summary of the global evolution previous to 
TP-AGB phase 

4.1. Some results 

Theoretical Hertzsprung-Russell diagrams (HRD) for the 
ensemble of our models are shown in Fig. [l], from the early 
PMS phase up to the beginning of the first thermal pulse 
(or He-shell flash) along the AGB phase. 

The lifetimes are given in Table |l| for the PMS, the H- 
burning phase (from the zero-age-main sequence to central 
H-exhaustion) , the RGB phase (from central H-exhaustion 
to central He-ignition) , the He-burning phase (up to cen- 
tral He-exhaustion), the early (E) AGB phase (up to the 
first thermal pulse) and the thermally pulsing (TP) AGB 
phase. This last duration includes our extrapolation pro- 
cedure up to the convective envelope exhaustion (see Sect. 
6.3). 

The occurrence of the first thermal pulse marks the end 
of the E-AGB phase and the beginning of the TP-AGB 
phase. This precise time, marked by the maximum surface 
luminosity just preceding the first major thermal instabil- 
ity, will be called to in the following. Table || presents some 
characteristics of our models at tg. 

4.2. Comparisons with other works 

For the main sequence duration, our results for Z = 0.02 
are in better agreement with Bressan ct al. (1993) predic- 
tions than with the slightly higher values given by Schaller 
et al. (1992). The central helium burning lifetime crucially 
depends on the amount of overshooting. As our stellar 
models are computed without overshooting, we can only 
compare our predictions at this phase with those of the 



Padova group calculated in the same conditions; we ob- 
tain very similar results. 

The other important quantity for the AGB phase is 
the core mass (i.e. the mass contained up to the HBS) 
at the end of the central He-burning phase. Indeed, its 
value determines the surface luminosity (through the core 
mass-luminosity relation; see e.g. Boothroyd & Sackmann 
1988a), the mass loss rate, and consequently the AGB 
phase duration. Last but not least, along the TP-AGB 
phase, the thermal pulses are stronger and the third 
dredge- up deeper when the core mass is higher. Our 3 Mq 
star with Z = 0.02 enter the AGB phase with a core mass 
of 0.542 Mq (see Table a value very close to that of 
the comparable model of Boothroyd & Sackmann (1988a) 
or Lattanzio (1986), also computed without overshooting. 
Always for a reference 3Mq star, the TP-AGB phase of 
the Straniero et al. (1995) models begins with a core mass 
of 0.53 Mq (2 % lower). At a comparable evolution stage 
however, the 5 Mq TP-AGB star of Vassifiadis & Wood 
(1993), also computed without any kind of extra- mixing, 
has a core mass ~ 4 % higher than ours. This could be 
due to the somewhat lower metallicity of their models 
{Z = 0.016). In conclusion, the core mass at the begin- 
ning of the TP-AGB phase, although very important for 
this phase, still depends on various uncertainties related 
to the treatment of convection during the previous evolu- 
tionary phases. 

5. Surface abundances prior to TP-AGB phase 

5.1. First and second dredge-up events 

During the first and second dredges-up, convective mix- 
ing and induced dilution modify the surface abundances 
which were unaltered until this phase (Iben 1964). Table 
H presents the resulting changes for the principal elements 
altered by H-burning. Let us emphasize the main points. 

A second dredge-up occurs in our > 4Mq models, also 
leading to surface abundance changes. For our > 5Mq 
models, the convective envelope penetrates so deep that 
it pushes down the H-He discontinuity. 

5.1.1. ^He 

■^He is produced at the beginning of the T Tauri phase 
from the initial ^H burning. On the main sequence, a peak 
of -^He appears, due to the competition of the reactions 
p (p , e+ i^e) ^H (p , 7) ^He and ^He (^He , 2p) ''He . During 
the first dredge-up, the surface abundance of "^He thus 
increases. However, the higher the stellar mass, the lower 
the '^He peak, and the lower the surface enhancement of 
■^He after the first dredge-up. 

5.1.2. CNO elements 

The first (and second) dredge-up phase (s) lead to a de- 
crease of the 12C/13C, 12C/14N, 15N/14N and i^O/i^O iso- 
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Table 1. Lifetimes in all the evolutionary phases for the modeled stars which undergo thermal pulses 









Z 


= 0.02 






Z = 0.005 






3M0 


4M0 


5M0 


6M0 


3M0 


4M0 


5M0 


At [PMS] (yr) 


1.20 lO'^ 


5.88 


10® 


3.70 10® 


2.19 10® 


1.10 lO'^ 


5.82 10® 


3.50 10® 


At [MS] (yr) 


2.80 10* 


1.34 


10** 


7.69 10^ 


5.09 10'^ 


2.5910* 


1.30 10* 


7.88 lO'^ 


At [RG] (yr) 


1.78 lO'^ 


6.39 


10® 


2.86 10® 


1.22 10® 


1.45 lO'^ 


5.1710® 


2.38 10® 


At [HB] (yr) 


1.55 10® 


5.26 


10^ 


2.64 lO'^ 


1.08 lO'^ 


9.3710^ 


4.30 10'^ 


2.59 lO'^ 


At [E-AGB] (yr) 


6.1610® 


2.88 


10® 


1.25 10® 


5.98 10^ 


4.29 10® 


1.69 10® 


6.41 10® 


At [TP-AGB] (yr) 


1.22 10® 


2.26 


10^ 


1.33 10^ 


9.12 10* 


9.03 10^ 


1.59 10^ 


7.97 10" 



Table 2. Main features concerning the internal and surface structure of our models at the end of the E-AGB phase. [He] {[H]) 
specifies values taken at the maximum energy production rate in the helium (hydrogen) burning shell. [H — He] refers to the 
intcr-shcll region (i.e. between the HeBS top and the HBS base) and [env] to the convective envelope bottom. Am is the mass 
thickness of a region 



Z = 0.02 Z = 0.005 





3M0 


4M0 


5M0 


6Mq 


3M0 


4M0 


5M0 


^0 (yr) 


4.739 10* 


2.034 10* 


1.123 10* 


6.624 10^ 


3.851 10* 


1.880 10* 


1.11710* 


L(L0) 


3921 


13274 


19367 


25743 


5420 


14761 


21479 


^(Rq) 


182 


398 


486 


554 


182 


339 


387 


AR[atmos.] (Rq) 


3 


15 


17 


18 


3 


10 


10 


M(Mq) 


2.92 


3.77 


4.70 


5.72 


2.81 


3.51 


4.31 


M(M0yr-") 


2.010"'' 


1.5 10"® 


2.5 10"® 


3.610"® 


3.610"^ 


1.8 10"® 


2.810"® 


Teff (K) 


3381 


3102 


3088 


3104 


3667 


3450 


3550 


loggeff (gcm~^) 


0.38 


-0.19 


-0.26 


-0.29 


0.36 


-0.08 


-0.10 


Lnucl/ L 


0.96 


0.87 


0.80 


0.71 


0.94 


0.79 


0.66 


Lhr/Lh 


0.09 


0.08 


0.07 


0.06 


0.06 


0.13 


0.10 


Vc 


20 


30 


37 


39 


20 


32 


39 


Tfnax (-^) 


2.11 10* 


3.48 10* 


3.97 10* 


4.60 10* 


2.44 10* 


3.91 10* 


4.54 10* 


rriTmax (Mq) 


0.320 


0.521 


0.643 


0.745 


0.373 


0.608 


0.736 


r[He] (Ro) 


0.00008 


0.00003 


0.00003 


0.00002 


0.00008 


0.00004 


0.00003 


m[He] (Mo) 


0.520 


0.732 


0.815 


0.894 


0.594 


0.814 


0.903 


Am[He] (Mq) 


0.08 


0.03 


0.02 


0.24 


0.04 


0.01 


0.23 


T[He] (K) 


1.43 10* 


1.61 10* 


1.64 10* 


1.68 10* 


1.37 10* 


1.6410* 


1.68 10* 


p[He] {gcm-^) 


13.8 


7.7 


5.9 


5.9 


11.2 


5.3 


4.6 


Am[H - He] (M®) 


4.0 10"^ 


4.710"" 


1.510"" 


6.010"® 


3.3 10"^ 


3.8 10"" 


1.2 10"" 


r[H] (Rq) 


0.00016 


0.00006 


0.00004 


0.00003 


0.00015 


0.00006 


0.00005 


m[H] (Mq) 


0.542 


0.738 


0.818 


0.895 


0.609 


0.818 


0.905 


Am[H] (Mo) 


1.210"^ 


3.2 10"" 


1.810"" 


5.510"^ 


9.710"" 


3.4 10"" 


2.210"" 


T[H] (K) 


5.0710^ 


6.43 10^ 


7.10 10^ 


7.77 lO'^ 


5.75 lO'^ 


7.00 10^ 


7.74 10^ 


p[H] (gcm-^) 


2.010"^ 


1.4 10"^ 


1.310"^ 


1.1 10"^ 


2.410"^ 


1.710"^ 


1.610"^ 


r[env] (Ro) 


0.0042 


0.0008 


0.0002 


0.0001 


0.0044 


0.0008 


0.0002 


m[enii] (Mq) 


0.544 


0.738 


0.818 


0.896 


0.612 


0.819 


0.905 


T[env] (K) 


2.68 10® 


5.72 10® 


1.58 10^ 


4.2710'' 


2.71 10® 


6.53 10® 


2.25 10'' 


p[enw] (gcm~^) 


1.1 10"^ 


4.1 10"^ 


5.7 10"^ 


0.9 


1.010"^ 


6.0 10"^ 


0.2 
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Fig. 1. Theoretical HRD, from the beginning of the PMS phase up to the beginning of the first thermal pulse along the AGB 
phase, for Z = 0.02 (left) and Z = 0.005 (right) 



topic ratios, the extent of which depends on both stellar 
mass and metallicity. 

— From an initial value equal to 83, the -'^^C/^'^C ratio 
decreases to ~ 18 — 19. The higher the stellar mass, 
the lower the final carbon isotopic ratio. 

This ratio slightly further decreases after the second 
dredge-up. 

— For Z ~ 0.02, the maximum depth of the convective 
envelope reaches the deepest step of the ^'^N profile, 
for > SMq stars. This implies that ^^C/^N, as well 
as ^^N/^^N have almost the same value in these mod- 
els after the first dredge-up. On the contrary, for the 
Z = 0.005 objects, the lower the stellar mass, the less 
efficient the dredge-up, and consequently, the less ^"^N 
is up-heaved to the surface. This explains the trend for 
the final ^^C/^'*N ratio versus stellar mass (see Char- 
bonnel 1994 for more details). 

After the second dredge-up, whatever Z , the convec- 
tive envelope of > 4Mq stars penetrates deeper for 
increasing stellar mass as attested by both ^^C/^^N 
and i^N/"N ratios. 

— The theoretical post dredge- up value of ^^O/^^O in- 
creases with the stellar mass. Moreover, since for Z = 
0.005 the deepening convective envelope of our 4 and 
5 Mq stars does not reach the regions where the ON 
cycle operated, the post dredge-up values of ^^O/^^O 
are higher than for Z ~ 0.02. 



After the second dredge-up however, ^^O enriched mat- 
ter is up-heaved in the convective envelope leading to 
very low ^Sq/^O. 
— The predicted post dredge- up value of ^^O/^^O shghtly 
decreases with increasing stellar mass. 



5.2. Comparisons with other works 

Our ^He surface abundance increases are in rather good 
agreement with the predictions by Bressan et al. (1993). 

Concerning the CNO elements, our results are in gen- 
eral good agreement with the predictions by Schaller et al. 
(1992), Bressan et al. (1993) and El Eid (1994). The only 
important discrepancy between theoretical predictions of 
the different groups for intermediate-mass stars concerns 
the post dredge-up ^^O/^^O ratio, which value highly de- 
pends on the adopted ^"^0 (p , a) ^^^N and "O (p , 7) ^^F 
reaction rates. We use the rates given by Landre et al. 
(1990), as do Schaller et al. (1992), and we obtain ^^Q/i^O 
ratios very similar to those of the Geneva group. On the 
other hand, Bressan et al. (1993) use the lower ^^O proton 
capture rates given by Caughlan & Fowler (1988), and ob- 
tain much lower -^^O/^^O ratios [see also El Eid (1994) for 
a discussion of the influence of the adopted ^^O destruc- 
tion rate on the resulting ^^O/^^O ratio]. 
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Table 3. Surlace abundances and isotopic ratios consecutive to the first and second dredge-up episodes (i.e. respectively at the 
top RG and time to defined in the text) 









Z 


= 0.02 






Z = 0.005 




3M0 


4Mq 


5Mq 


6M0 


3M0 


4M0 


5Mq 


^He 


initial 


3.22 10"^ 


3.1910"^ 


3.1310"^ 


3.1210"^ 


2.71 10"^ 


2.71 10"^ 


2.66 10"^ 




top RG 


1.55 10"" 


1.05 10"* 


8.52 lO""^ 


7.32 10"^ 


1.71 10"" 


1.28 10"" 


1.08 10"" 




at to 


1.54 10"" 


1.03 10"" 


7.79 lO""^ 


6.65 10"^ 


1.66 10"" 


1.13 10"" 


8.56 lO"'""' 


*He 


initial 


0.2932 


0.2932 


0.2932 


0.2932 


0.2495 


0.2495 


0.2495 




top RG 


0.3047 


0.3039 


0.2995 


0.3017 


0.2585 


0.2505 


0.2498 




at to 


0.3054 


0.3052 


0.3367 


0.3421 


0.2624 


0.2629 


0.3014 




initial 


9.90 10"^ 


9.90 10"^ 


9.90 10"^ 


9.90 10"^ 


2.4710"^ 


2.4710"^ 


2.47 10"^ 




top RG 


1.4410"^" 


1.44 10"^° 


1.43 10""^ 


1.44 10"^° 


3.1410"" 


3.51 10"" 


3.86 10"" 




at to 


1.20 10"^" 


9.85 10"" 


7.39 10"" 


7.53 10"" 


1.49 10"" 


1.68 10""^ 


3.59 10"^** 




initial 


83.1 


83.1 


83.1 


83.1 


83.1 


83.1 


83.1 




top RG 


19.3 


19.0 


18.8 


18.5 


18.7 


18.4 


18.3 




at to 


19.1 


18.6 


17.8 


17.5 


18.1 


17.5 


17.4 


12C/14N 


initial 


2.745 


2.745 


2.745 


2.745 


2.745 


2.745 


2.745 




top RG 


0.725 


0.714 


0.769 


0.735 


0.636 


0.888 


1.100 




at to 


0.712 


0.684 


0.565 


0.538 


0.581 


0.610 


0.496 




initial 


253 


253 


253 


253 


253 


253 


253 




top RG 


1243 


1281 


1210 


1279 


1397 


1033 


856 




at to 


1270 


1346 


1654 


1774 


1546 


1510 


1841 


12(./16o 


initial 


0.316 


0.316 


0.316 


0.316 


0.316 


0.316 


0.316 




top RG 


0.205 


0.207 


0.210 


0.211 


0.194 


0.208 


0.229 




at to 


0.204 


0.205 


0.206 


0.207 


0.191 


0.192 


0.197 


160/170 


initial 


2467 


2467 


2467 


2467 


2467 


2467 


2467 




top RG 


314 


409 


562 


617 


258 


869 


1814 




at to 


311 


391 


508 


579 


245 


323 


358 


I60/I80 


initial 


442 


442 


442 


442 


442 


442 


442 




top RG 


589 


582 


577 


577 


616 


573 


522 




at to 


591 


588 


588 


588 


623 


619 


611 



5.3. Comparisons with observations 

In our intermediate-mass stars, the theoretical post 
dredge-up values of the carbon isotopic ratio are slightly 
lower than the observations by Gilroy (1989) in galactic 
cluster giants. 

Red giants in this evolutionary phase present ^^O/^^O 
ratios between 300 to 1000 and ^^O/^^O ratios in the range 
400 to 600 (Harris et al. 1988; Smith & Lambert 1990a). 
This is in very good agreement with our predictions. In 
particular, our prediction is in perfect agreement with the 
observed value in a UMa (Harris et al. 1988), which es- 
timated mass is roughly 3 Mq . This point is in favor of 
high ^^O destruction rates. 

This agreement between standard models and obser- 
vations confirms that no extra-mixing (diffusion, rotation- 
induced mixing, . . . ) is expected to occur in intermediate- 
mass stars, as already discussed in Charbonnel (1994, 
1995). 



5.4. Internet Tables 

Very complete tables containing informations about the 
center and surface evolution of the structure and chemical 
composition for our 3, 4, 5, 6 and 7Mq models with Z = 
0.005 and 0.02, from the beginning of the PMS phase up 
to the end of the E-AGB phase, can be found through the 
internet network at the following address and path: 

pttp: / /gag.observ-gr.fr/liens/starevol/evol.htm| 

More specifically, we present, for each star, the evolu- 
tion of 

— the age, mass ,mass loss rate, surface luminosity, bolo- 
metric magnitude, effective temperature, total radius, 
surface density and gravity; 

— the central temperature, density, pressure, degeneracy 
parameter and the mass fraction of ^H, ''He, ^^C and 
160; 
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— the mass and radius positions of the central convective 
core boundary and of the bottom of the convective 
envelope as well as its temperature and density; 

— the nuclear luminosities of the H- and He-burning 
shells; 

— the mass and radius positions of the H- and He-burning 
shells; 

— the surface mass fractions of ^H, "^He, ^He, ^Li, ^Be, 
12c, 13c, i^N, 15N, 160, 1^0, iSQ, i^F, 2"Ne, ^iNe, 22Ne, 

23Na, 24Mg, 25Mg, 26Mg, 26A1S, 27A1 

6. structural properties of the thermal pulses 

6.1. Evolutionary properties of the thermal pulses as a 

function of M and Z 

As the stars are ascending the AGB, the features of the 
thermal pulses evolve, and some of the main characteriz- 
ing quantities of these events reach "asymptotic values" . 
The most important feature is the thermal pulse inten- 
sity that can be defined in terms of the HeBS luminos- 
ity i^g"^ reached at the maximum of the He-shell flash, 
^maa; jg q^j^g Jqw during the first thermal pulses, but it 
rapidly increases with core mass (or flash number). The 
increasing rate is greater for lower mctallicitics or higher 
stellar masses. After typically 10 (5) thermal pulses in our 
Z = 0.02 (0.005) models, the thermal pulse intensity is 
very high and only slightly increases from pulse to pulse, 
monotonously with the core mass. This second part of the 
TP- AGB phase is usually called the "asymptotic regime" 
or "full amplitude" of the thermal pulses. 

For each of the modeled stars, we have computed at 
least four full amplitude thermal pulses. Let us briefly dis- 
cuss the main evolutionary properties of these thermal 
runaways. 

The thermal runaways are stronger inside stars of lower 
Z, which facilitates dredge-up events (sec Sect. 6.2). The 
intensities of our He-shell flashes are comparable to other 
published values for stars of about the same core masses 
(see e.g. Boothroyd & Sackmann 1988b and Table 4). 

While the stars evolve along the TP- AGB, the HeBS 
and HBS advance in mass, and the mean inter- shell mass 
slightly decreases with time. From one pulse to another, 
the top of the convective tongue at its maximum extension 
gets closer from the HBS, even if contact between both 
regions has never been found due to the fact that the HBS 
is still active at that time. Moreover, for all the thermal 
runaways we have computed, the overlap in mass between 
two successive thermal pulses is ranged between 0.6 and 
0.4 for AGB stars of increasing total mass. This overlap is 
slightly larger for our lower Z models. 

The maximum (minimum) temperature T™°^ ^jmin^ 
and density p™"^ (p™*") reached at the base (top) of a 
thermal pulse both increase (decrease) from pulse to pulse. 
For a given He-shell flash, T™"^ increases with stellar mass 
and metallicity, while p™"^ decreases with stellar mass but 



increases with metallicity. In the asymptotic regime, the 
growth rate of T™"^ also tends to an asymptotic value. 
By comparison with other works for the 3 Mq (Z = 0.02) 
star, we found that our T™"^ are typically 4 % higher than 
those of Boothroyd & Sackmann (1988b) and Straniero et 
al. (1996). This could be mainly due to (i) different pre- 
scriptions to compute Sgrav, (a) different envelope masses 
due to different mass loss rates and (Hi) different MLT 
treatments (see Sect. 2.1.4). The temperature and density 
differences between the base and top of each convective 
tongue are rather large. This point is important for nucle- 
osynthesis purposes (see Sect. 7). 

The duration of the convective tongue associated with 
each thermal pulse (Aip) and, to a lower extent, that of 
the inter-pulse phase (At^p), are crucial quantities for the 
nucleosynthesis too. The stronger a He-flash, the shorter 
it is and the quickest it deactivates. Both quantities are 
thus smaller for higher core masses and lower metallicities. 
Moreover, Atp tends to decrease as the stars evolve along 
the TP- AGB. Afjp also decreases after the first He-shell 
flashes, but it then reaches an almost constant value dur- 
ing the asymptotic regime, only very slightly decreasing 
with the increasing core mass. 

6.2. Third dredge-up 

Along the TP-AGB phase, the bottom of the convective 
envelope can penetrate deep inside regions that have been 
nuclearly processed just after a thermal pulse, if it is 
strong enough (in order to give rise to a large enough en- 
tropy increase in the inter-shell region). Contrarily to the 
flrst and second dredges-up, this so-called third dredge- up, 
when it occurs, successively mixes up to the surface (i) ma- 
terial that has experienced H-burning in the thin HBS but 
also (ii) part of the region where the thermal pulse nucle- 
osynthesis operated just before. Both regions give rise to 
very different and sometimes opposite chemical pollu- 
tions of the convective envelope (see Sects. 7.1 and 7.6). 

In our Z = 0.02 models, the HBS region is already 
dredged-up from the 14th, 5th and 3rd thermal pulse 
for the 3, 4 and 5M0 stars, respectively, while for the 
Z = 0.005 models, such a mixing already occurs from the 
2nd and 1st thermal pulse for the 3 and 4Mq stars, re- 
spectively. The 6M0 (5M0) star with Z = 0.02 (0.005) 
continuously mixes H-burning products up to the surface 
as it cxpcric;nc;c;s H-burning at the bottom of its convective 
envelope (the so-called "Hot-Bottom Burning" or HBB; 
see Sect. 7.1). This kind of dredge-up already modifies the 
surface composition for some species (like ^^Al^; see Sects. 
7.4 and 7.6). 

The complete third dredge-up (one usually simply calls 
the third dredge-up, 3DUP hereafter), i.e. reaching the 
upper part of the previous thermal pulse region, is more 
difficult to obtain numerically. Following models (Wood 
1981; Iben 1983; Lattanzio 1987; Boothroyd & Sackmann 
1988c), the dredge-up of thermal pulse material is eas- 
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ier for stars of lower metallicity Z . This is supported by 
observations (see e.g. Blanco et al. 1978). In our mod- 
els, we just found a 3DUP from the 9th thermal pulse of 
the 4Mq with Z ~ 0.005 star. However, based on the 
evolution of the post-thermal pulse entropy profile from 
pulse to pulse, we can predict that such a dredge-up will 
also occur, within the next two thermal pulses, for our 
5 and 6Mq stars with Z = 0.02 and for the 5Mq star 
with Z = 0.005 (see Table ^ for the number of computed 
thermal pulses for each object). Some additional thermal 
pulses are still needed for lower mass objects. Several years 
ago, the 3DUP was found to occur too late along the TP- 
AGB phase for AGB stars of decreasing mass and/or in- 
creasing metallicity (see e.g. Iben 1976, Fujimoto & Sug- 
imoto 1979). This was a common problem (see Lattanzio 
1989 for a review). Recent works showed that with the new 
OPAL opacities, the 3DUP can now be obtained even for 
low-mass Pop I objects (see e.g. Frost & Lattanzio 1996 
or Straniero et al. 1995). 

One of the clearest signature of the 3DUP is the signifi- 
cant ^^C pollution of the convective envelope, brought into 
the inter-shell region by the convective tongue of a ther- 
mal pulse (see Sect. 8.1). As the thermal pulse intensity 
(and consequently the probability of 3DUP occurrence) 
increases with time, as do the core mass and the total lu- 
minosity, the models tend predict the formation of carbon 
(C) stars (for which the surface ^^C/^^0 becomes greater 
than unity thanks to the repetitive 3DUP events) at too 
high luminosities, i.e. to late during the TP-AGB phase, 
especially for solar metallicity stars. Solutions to increase 
the depth of the convective envelope during a dredge-up, 
e.g. by treating in some way semi-convection (see e.g. Lat- 
tanzio 1986) or by invoking undershooting (see Alongi et 
al. 1991) have been suggested. They indeed help the 3DUP 
to occur but do not yet solve completely the problem. Fun- 
damentally, this could be due to our bad knowledge of con- 
vection inside stars. Very recently, Straniero et al. (1996) 
suggested that by considerably increasing the number of 
mass shells and especially the number of time steps to 
model thermal pulses, the 3DUP occurs earlier during the 
TP-AGB phase. This consists in a real numerical progress 
in that field. The fact we find 3DUP latter than them 
along the AGB phase can be due to our lower number of 
time steps (typically by a factor of 3). 

6.3. Towards the end of TP-AGB phase 

6.3.1. Core Mass-Luminosity relation 

It was first discovered by Paczynski (1975) that during 
the asymptotic TP-AGB phase, there is a linear rela- 
tion between the surface luminosity L and the core mass 
Mc = m[H] (see Table H) which is only slightly dependent 
on the total mass of the AGB star. Blocker & Schonberner 
(1991) however found that this kind of relation does not 
exist for the most massive AGB stars for which the bot- 



tom of the convective envelope goes down to the HBS so 
that part of the H-burning occurs inside the convective 
envelope (HBB). 

^From our models for the last four thermal pulses of 
the 3 and 4M0 with Z ^ 0.02 and the 3M0 with Z = 
0.005, we have derived a core mass-luminosity relation in 
the form 



L 



M 
/ M 



0.02 

Z 



7 r 



Mc 







(5) 



where a = 1233L0, (3 = 3.566, 7 = 0.519, 5 = 0.1057, 
£ = 1.312 and C = 0.203. It is in good agreement with 
other Mc - L relations, for 0.50 < Mc/Mq < 1.10 and 
0.005 < Z < 0.02 (see e.g. Lattanzio 1986 or Boothroyd 
& Sackmann 1988a). The more massive stars experience 
HBB and do not fit, consequently, to the same relation. 
In Table ^, we however give a linear relation of the form 
L/Lq = a{Mc/MQ-(j) for the 5M0 (4M0) AGB star 
with Z = 0.02 (0.005) that do not experience a very strong 
HBB. Let us stress that during the TP-AGB phase, L suf- 
fers very important variations during and after each ther- 
mal pulse. Frequently, the core mass-luminosity relation 
is derived by taking solely account of the maximum sur- 
face luminosity preceding each thermal pulse. However, 
the probability to observe an AGB star just at that time 
is very weak (~ Atp/AUp < 1%; see Table ^. We thus 
preferred to derive our relations with a more realistic sur- 
face luminosity that is averaged over each pulse-inter-pulse 
cycle. 

On the other hand, it is well known that the rate of 
mass advance of a burning shell (here Mc for the HBS) is 
related to the rate of nuclear energy production by 



dMc 
dt 



L 



H 



XEh 



(6) 



with Lh being the luminosity produced by the HBS, X 
the mean hydrogen mass fraction over this burning shell 
and Eh = 6.35 10^^ erg g~^ the amount of energy pro- 
duced by the conversion of 1 g of -'^H to ''He through the 
CNO bi-cycle. Along the TP-AGB, most of the surface 
luminosity L is produced by the HBS. As indicated in 
Table |^, the HeBS contributes much less than the HBS 
to the nuclear luminosity (except during a thermal run- 
away). Gravitation also constitutes a relatively small en- 
ergy source, especially for the less massive AGB stars. In 
conclusion, we can write Lu — £,L with 0.6 < ^ < 0.9 de- 
pending on M. By eliminating Lh in Eq. (6) with the core 
mass-luminosity relation, one finds that the core mass Mc 
increases exponentially with time during the asymptotic 
TP-AGB phase. The time scale for the core mass growth 
is - 10^ (- 10^) yr for a 3M0 (6M0) AGB star, that 
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is, interestingly, comparable to the total duration of the 
TP-AGB phase. Table ^ demonstrates that Eq. (6) pre- 
dicts rates of HBS advance in mass that are very similar 
to those coming from our numerical models, except for the 
6M0 (5M0) object with Z = 0.02 (0.005), probably due 
to the strong HBB. 

6.3.2. Extrapolation procedure 

One of our ambitions is to extrapolate our predictions 
coming from the complete evolution models not only for 
the global surface properties but also for the chemical com- 
position. Eqs. (5) and (6) already allow to extrapolate 
rather accurately the core mass and surface luminosity 
from our last computed models towards the end of the 
TP-AGB phase. For the 6Mq (SM©) star with Z = 0.02 
(0.005), independent fits of L and Mc as a function of 
time have been calculated. We also need fits for the sur- 
face radius [in order to follow the mass loss rate increase 
through Eq. (1)] as well as the temperature and density 
at the base of the convective envelope (to follow the HBB 
nucleosynthesis). We found, on grounds of the last four 
computed thermal pulses for each star, very good fitting 
relations for these quantities, which are presented in Ta- 
ble ^. Note that they are linear with time during the last 
part of the TP-AGB phase, contrarily to Mc and L that 
increase exponentially during this phase. 

The main interest of such an extrapolation procedure 
is that very detailed stellar evolution computations (in- 
cluding a very complete description of the nucleosynthe- 
sis) would take prohibitive computation time. However, we 
had to make two important approximations (compared to 
full evolution models) that must be mentioned. 

— The successive 3DUP events progressively pollute the 
convective envelope in ''He, '^C and heavier elements 
to the detriment of ^H. This increases the mean molec- 
ular weight (or Z) and the mean opacity in the sub- 
photospheric region, probably leading to an overesti- 
mation of the surface luminosity. 

— Very few computations are existing for the very end of 
the TP-AGB phase (Vassihadis & Wood 1993; Blocker 
1995), that cannot maintain to give accurate predic- 
tions concerning the evolution of the surface isotopic 
composition. One of their interests is to indicate that 
the last few thermal pulses seem to have different prop- 
erties than the full amplitude ones (e.g. longer inter- 
pulse durations). On the other hand, it has been noted 
by Boothroyd & Sackmann (1988a) that the increasing 
mass loss rate and Z possibly prevent the 3DUP to oc- 
cur after these last thermal pulses. As a consequence, 
assuming the full amplitude regime up to the AGB tip 
probably makes our following predictions quite indefi- 
nite concerning the very end of the TP-AGB phase. 

Such extrapolation procedures have already been re- 
alized by Groenewegen & de Jong (1993 and following 



papers). These synthetic evolution models are very pow- 
erful to compare observations with computed models con- 
cerning the evolution of the global surface properties of 
AGB stars. However, they contain shortcomings that do 
not allow to predict in detail the evolution of the surface 
isotopic composition. Mainly, (i) only the most abundant 
elements are followed (''He, '^C, '^C, '''N and '^O) and (n) 
the HBB nucleosynthesis is not followed in detail whereas 
it becomes a very influent process concerning the convec- 
tive envelope composition towards the TP-AGB tip. Re- 
cent progresses have been made by Busso et al. (1995) and 
Lambert et al. (1995) to follow in detail the nucleosynthe- 
sis as late as possible along the AGB phase. 

In Sect. 8, we both present the nucleosynthesis encoun- 
tered in our full evolution models and the extension of the 
predictions concerning the convective envelope composi- 
tion towards the end of the AGB phase corresponding to 
the seven computed stars. For this extrapolation, we pro- 
ceed as follows. 

1. From the last computed inter-pulse phase for each star, 
we integrate, during the successive inter-pulse phases, 
the evolution of the nucleosynthesis occuring inside the 
convective envelope by solving Eqs. (4), suited for the 
HBB, for each nuclide while imposing a temperature 
and density increases at the bottom of the convective 
envelope given by the fits presented in Table ^. We 
took the inter-pulse durations of Table indeed, they 
are only very slightly decreasing with the core mass 
increase along the asymptotic TP-AGB phase. 

2. Since the 3DUP is almost ready to occur in our models 
(see Sect. 6.2), we mix the envelope with the HBS and 
inter-shell regions recurrently, every Ai^p. To quantify 
the adopted depth of each 3DUP, let us remember the 
definition of the usual parameter 



where IS.Mdijp represents the amount of mass that 
is dredged-up (difference between the normal depth 
of the convective envelope and its maximum down- 
wards penetration during a dredge-up) and AMc the 
increase in mass of the HBS between two successive 
thermal pulses. We adopted a 3DUP depth correspond- 
ing to A = 0.6 for all the stars. This value has been 
suggested by Groenewegen & de Jong (1993) because 
it reproduces rather well the carbon star luminosity 
function of the LMC (that cannot be fitted on grounds 
of self-consistent evolution models only; see Sect. 6.2). 
During such a dredge-up, we assumed that the mate- 
rial coming from the convective tongue of a thermal 
pulse has abundance profiles corresponding to end of 
the last computed thermal pulse. Of course, we cannot 
do better without further computing thermal pulses 
with detailed nucleosynthesis up to the end of the TP- 
AGB phase that is, again, not yet realistic. However, 
this appears to be rather justified for many important 
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Table 4. Some features of our TP- AGB models at the end of the complete stellar evolution computations. We indicate the total 
number Hp of computed thermal pulses. The core mass (Adc) - luminosity (L) relations L — a{Adc — P) indicated for the 5 Mq 
(4M0) star with Z = 0.02 (0.005) is marginally valid as these stars already experience significant hot-bottom burning. No such 
a relation is given for the 6Mq (5Mq) object with Z = 0.02 (0.005). We refer to the text for the fitting core mass-luminosity 
relation found for the other stars. M^nv is the envelope mass. We also indicate the intensity of the last computed thermal 
pulse (in term of the maximum luminosity L^e^ produced by the HeBS during this thermal runaway), the duration of the 
last intor-pulsc (Atip) and pulse (Atp), the time Atp^ elapsed between the maximum extent of the convcctive tongue of the 
last thermal pulse and its disappearance, the maximum mass extension of the convective tongue (Amp) and the mass overlap 
between the last two thermal pulses {Ampr). We finally mention the occurrence of the third dredge- up 
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Table 5. Comparison between predictions for the rate of HBS advance in mass from Eq. (6) and those coming from our 
computed models. The mean hydrogen mass fraction X and HBS luminosity Lh are taken during the inter-pulse phase of the 
last computed thermal pulse for each star. dMc jdt is correspondingly calculated over the same inter-pulse phase 
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chemical elements as they reach almost constant abun- 
dances after the few last computed thermal pulses (see 
discussion in Sect. 7). 

Taking a constant value for A is probably not strictly 
right. As far as the surface composition is concerned, in- 
creasing A more efficiently pollutes the convective envelope 
from pulse to pulse. 

There is another quantity that greatly influences the 
evolution of the envelope composition: the mass loss rate. 



Indeed, decreasing the mass loss rate allows to increase 
the total number of thermal pulses (and consequently the 
total number of 3DUP episodes) before the end of the TP- 
AGB phase. We had to increase our r] values in Eq. (1) 
for each star, in order to reach the highest mass loss rates 
(called super-winds) observed among very evolved AGB 
stars (see the discussion in Sect. 2.1.5: sec also Zuckerman 
et al. 1986; Wannier & Sahai 1986; Knapp et al. 1989): we 
adopted ?? = 5, 10, 15 and 20 for the 3, 4, 5 and 6 Mq AGB 
models, respectively, whatever Z. As the exact way follow- 
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Table 6. Fitting parameters from the last four computed asymptotic thermal pulses towards the end of the TP-AGB phase, 
for the mean total radius and bottom temperature and density of the convective envelope. For each quantity A, the relation 
reads A{t) — a{t — to) + 13, with time in yr 
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ing which the mass loss rate increases with time along the 
asymptotic TP-AGB phase is still questionable by more 
than a factor of two (see Sect. 2.1.5), the resulting uncer- 
tainty on the evolution of the surface composition up to 
the AGB tip dominates the effect of A (see our discussion 
in Sect. 8.1, where our mass loss rates are varied by a 
factor of four). 

In Table 0, we present the global mass loss increases we 
obtained, for the seven TP -AGB stars, from the beginning 
of the TP-AGB phase up to its end (i.e. including our 
extrapolation computations). 

6.3.3. Initial mass-Final mass relation 

One of the most important constraint for evolutionary 
models towards the AGB tip is the reproduction of the 
observed correlation between the initial mass of stars (Mi) 
and their final one (M/), i.e. their core mass (become a 
hot white dwarf) when the convective envelope has been 
completely removed after the planetary nebula (PN) ejec- 
tion. We present in Fig. ^ the comparison between ob- 
servations analyzed by Weidemann & Koester (1983) and 
Weidemann (1987) and the final core masses we obtained 
after our evolutionary plus extrapolations computations. 
We note that agreement with observations is slightly bet- 
ter than for Vassiliadis & Wood (1993), especially around 
Mi = 3 Mq where they only reproduce the upper lim- 
its. On the contrary, the final core mass obtained by 
Boothroyd & Sackmann (1988b) is rather low compared to 
observations. Note we have not taken into account of the 
core mass reduction following each 3DUP event. However, 
as explained below, such a reduction, even if it increases 
the agreement with the observed relation, is rather neg- 
ligible compared to other uncertainties (core mass at the 
first thermal pulse, thermal pulse number depending on 
the adopted mass loss rate, . . . ). 

We observe that changing the mass loss rate by a factor 
of four during the asymptotic TP-AGB phase (see also 
Sect. 8.1 concerning surface abundance predictions) does 
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Fig. 2. Relation between the initial mass of stars and the fi- 
nal white dwarf (WD) mass at the end of the AGB phase. 
Observational final masses are derived either from the surface 
gravity of the WD (filled circle) or from their radius (open 
circles); lines are connecting both determinations for identical 
objects. Our predictions are the black and gray rectangles, for 
the Z — 0.02 and Z = 0.005 modeled stars, respectively. Their 
extension is the result of varying the mass loss rate by a factor 
of four during the asymptotic TP-AGB phase 
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Table 7. Increase of the mass loss rates during the whole TP-AGB phase 
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not significantly change the agreement with observations. 
We claim that this is due to the fact that while almost 
50% of the C-0 core mass is built during the central He- 
burning phase, roughly 35 to 45% (in mass) are added 
during the E-AGB phase. More precisely, during the whole 
TP-AGB phase, we found that our core masses increased 
by 14 (12), 7 (4), 5 (3) and 4% for the 3, 4, 5, and 6Mq 
stars with Z — 0.02 (0.005), respectively. The final core 
mass is consequently almost decided before the super-wind 
occurs and this observed relation is not very compulsive 
for the empirical mass loss rate to adopt during the TP- 
AGB phase. 

7. Nucleosynthesis during TP-AGB phase 

^From the pioneer work of Iben (1975), one realized that 
inside TP-AGB stars, a very rich nucleosynthesis occurs. 
During the recurrent third dredge-up events, the syn- 
thesized nuclides are mixed inside the convective enve- 
lope and then ejected into the interstellar medium (ISM) 
through the strong mass loss. In conclusion, TP-AGB 
stars play a crucial role for the chemical evolution of 
galaxies, especially for some elements like "^He, ^Li, '^^F, 
^^Al, nuclides heavier than ^^Fe, ... It is the reason why 
we will give, in Sect. 8.3, the yields resulting from our 
intermediate-mass star computations. 

7.1. Overview of the various nucleosynthesis sites 

Let us analyze the various nucleosynthesis processes oc- 
curing inside intermediate-mass TP-AGB stars. In fact, 
three very different sites have to be distinguished. 

7.1.1. Burning shells 

As indicated in Table |^, most of the nuclear energy release 
on the AGB is due to the HBS, mainly through the CNO 
bi-cycle. More precisely, the six most important nuclear 
reactions that energetically supply the star are ^^C (p , 7) 
(- 26%), i4N(p,7) (~ 24%), i5N(p,a) (~ 17%), 
(/?+) (- 9%), (/3+) (^ 7%) and ^^C (p , 7) (- 6%), ah 
occuring inside the HBS. These numbers are quite similar 
(within a few %) for all our modeled stars. Whatever the 
AGB star total mass, the HBS slightly depletes ^^C and 
i^N and significantly depletes ^^O, ^^O and ^^F, while i**N 



is produced. The case of ^'^C is more difficult: while its 
abundance is very slightly increased in the HBS of 3 Mq 
stars (whatever Z), it is somewhat decreased in the HBS 
of 6 Mq stars. 

The rather high temperature inside the HBS also al- 
lows the activation of the NeNa and MgAl chains, even 
if they are much less influent for the global energetics. 
The most interesting consequences of the NeNa chain are 
the production of ^^Ne and ^'^Na in < 5Mq stars and 
their slight destruction in more massive AGB stars. The 
MgAl chain is very sensitive to the AGB total mass. In 
< 5 Mq objects, ^^Mg is destroyed to the benefit of ^^Al 
(roughly half in ^^Al^ and half in ^6^1™ that quite instan- 
taneously decays in ^^Mg) ; a small quantity of ^^Mg is also 
destroyed in ^^Al. In more massive AGB stars, the HBS 
temperature is high enough to substantially produce ^^Mg 
through the ^^Mg (p , 7) reaction. This in turn leads to a 
much important ^^Al^ production. However, as the HBS 
mean temperature always increases with time, ^^Al^ be- 
gins to be partially destroyed by (p, 7) reactions late along 
the TP-AGB. We refer to Forestini et al. (1991) for the 
first detailed study of ^^Al^ production inside AGB stars, 
by radiative proton burning (i.e. in the HBS). In Guelin 
et al. (1995), more details can be found, including the 
possibility to synthesize ^^Al^ by proton burning at the 
bottom of the convective envelope (i.e. by HBB). It was 
however a parametric study. First predictions concerning 
the enrichment of the AGB star surfaces, on grounds of 
full evolutionary models are presented in Sects. 7.6 and 
8.1. 

The comparatively small amount of nuclear energy 
supplied by the HeBS comes from the 3a (~ 70%) and 
^^C {a , 7) (~ 29%) reactions. Consequently, the HeBS 
mainly produces ^^C and, to a lower extent, "'^^O, while 
during the previous central He-burning stage, the ^^O pro- 
duction was preponderant. At the beginning of a thermal 
pulse, it is the 3a reaction that becomes the dominant 
energy source (see also Sect. 7.1.3 for the thermal pulse 
energetics) . 

During the long inter-pulse phase that follows a 3DUP 
episode, the ^^C coming out of the HBS can be burnt in the 
inter-shell region, close to the HeBS top, mainly through 
the ^^C {a , n) ^^O reaction. This reaction is highly ener- 
getic. The first published work reporting the existence of 
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this inter-shell ^'^C burning is Straniero et al. (1995). The 
produced neutrons are then captured by light species up 
to ^^Fe mainly, but part on them are also used to syn- 
thesize s elements. By decreasing order of amplitude, the 
nuclear energy liberated in the inter-shell region during 
the inter-pulse phase comes from these (n, 7) reactions 
up to 5^Fe (~ 29%), "c(a,n) (- 21%), the so-called 
s- process reactions (- 12%), ^SAls (n , p) ^SMg (~ 7%) 
and ^''N (n , p) ^''C 2.5%). Note the important nuclear 
energy released by the neutron captures. This energy pro- 
duction creates a third burning shell between the HeBS 
and HBS, even if its nuclear activity is considerably lower 
than that of both the HeBS and HBS. It continues to ex- 
ist until ^^C has been destroyed there. This happens, in 
general, almost when the following thermal pulse begins. 
As a consequence, when the convective tongue of a ther- 
mal pulse engulfs the inter-shell matter, the amount of ^^C 
that is mixed is reduced compared to what was previously 
thought, i.e. Am[H-He]X{^^C) [H-He], where [H-He] 
designates, as before, the inter-shell region. We show, in 
Fig. ^ the nuclear energy production rate {Snuc) profile in 
the case of our 3 Mq AGB star. This clearly demonstrates 
the three-burning shell structure. Such a structure is also 
found in the more massive AGB stars. Let us emphasize 
that this third burning shell could at least partially ex- 
plain the production of the s nuclides (see Sect. 7.5) in 
a radiative zone. This has been suggested by Straniero et 
al. (1995), even if their amount of ^'^C in the inter-shell 
was artificially increased in order to study the s— process 
nucleosynthesis. 

7.1.2. Hot-Bottom Burning 

If the temperature T^b at the base of the deep convective 
envelope is high enough (i.e. > 20 10^ K), i.e. if the base of 
the convective envelope reaches the HBS, proton burning 
occurs inside it. This HBB occurs in models with initial 
masses > 5Mq for Z = 0.02 and > 4M0 for Z = 0.005. 
These lower mass limits and this Z dependence are com- 
mon features of all the AGB models (see e.g. Boothroyd et 
al. 1993). As mixing up to the cool surface is very efficient 
(the mean turn-over time tco nv of the convective enve- 
lope is roughly 0.5 yr), the effective nuclear reaction rates 
are lower than in radiative H-Burning, leading to specific 
nucleosynthesis signatures. 

For Tcb > 40 10<^ K already, ^Li can be produced from 
''He following a scenario that cannot operate in radiative 
H-burning (see Sect. 7.2 for more details). This Li synthe- 
sis has been qualitatively suggested by Cameron & Fowler 
(1971) and quantitatively first demonstrated on grounds 
of evolutionary models by Sackmann & Boothroyd (1992). 

For Tcb > 60 10^ K, the CN cycle begins to operate par- 
tially, i.e. "'^^C is converted to ^^C that accumulates. As first 
investigated by Boothroyd et al. (1993), this conversion 
can efficiently reduce the ratio and eventually, if 
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Fig. 3. Profile of the nuclear energy production rate per unit 
mass, Enuc, in the region of the burning shells, in the case of the 
3Mq {Z — 0.02) star, during the inter-pulse phase between 
the 18th and the 19th thermal runaways. The dotted curve 
corresponds to a time just before the 19th thermal pulse, while 
the solid one is plotted roughly at the middle of the inter-pulse 
phase (see the text for a discussion) 

HBB is strong enough, this can prevent the formation of 
C stars (see also Sects. 8.1 and 8.2). 

If Tcb > 80 10^ K typically, the CN cycle almost op- 
erates at equilibrium through the whole envelope, leading 
to a ^^C/^'^C isotopic ratio lower than 10 (depending on 
Tcb)- In our models, this happens for the 6Mq (5Mq) 
AGB star with Z = 0.02 (0.005) and also, later along the 
TP-AGB, for our 5 Mq star with Z = 0.02 (see Sect. 8.1). 
At very high Tcb, other reactions begin to occur; in partic- 
ular, in extreme cases (i.e. the end of the TP-AGB phase 
of our 6M0 models with Z = 0.02), '^^AV can eventu- 
ally be directly produced inside the convective envelope 
through the ^^Mg (p , 7) reaction on ^^Mg supplied by 
24Mg(p,7) (see Sect. 8.1). 

7.1.3. Convective tongue of thermal pulses 

The thermal runaway, inside the HeBS, leading to form a 
convective tongue is principally due to the 3a reaction ex- 
citement. When this convective tongue finally penetrates 
the inter-shell region, the HBS ashes are ingested, i.e. not 
only ^He but also ^^C, ^''N and ^^Al^ mainly. As these ele- 
ments are partially burnt in the bottom of the inter-shell 
region during the inter-pulse phase (globally, we find that 
half of the ^■'C is so destroyed), they mainly survive in 
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the upper part of the inter-sheU. They are consequently 
engulfed by the convective tongue when it reaches its max- 
imum extension. At that time, the 3a reaction is already 
dominant for the energetics. We consequently find that 
the energetics associated with the ^'^C burning through 
the ^'^C (a , n) reaction as well as the neutron captures 
that follow is rather negligible (less than 1 %). This is 
also found by Straniero et al. (1995). This conclusion could 
however be modified if more ^^C was synthesized in the 
inter-shell region, as required for the s— process (see Sect. 
7.5). The important energetic effects reported by Bazan & 
Lattanzio (1993) concerning the ^'^C burning inside ther- 
mal pulses was probably due to (i) the fact they did not 
give account of partial inter-shell ^^C burning and (ii) the 
very high ^"^C amount artificially ingested into their ther- 
mal pulses in order to study the consequent s— process. 

Part of produced neutrons participate to (n, p) reac- 
tions, mainly on ^®AF and, to a lower extent, on ^''N. 
The resulting protons, at these high temperatures char- 
acterizing He-burning during a thermal pulse (> 200 10^ 
K at the base of the convective tongue), are rapidly cap- 
tured by various elements (see Sects. 7.3 to 7.5 for more 
details). Let us stress that the simultaneous presence, in- 
side a mixed He-burning region (the convective tongue) , of 
neutrons and protons leads to a very specific (and unique) 
nucleosynthesis. 

We will analyze in more details the nucleosynthesis 
processes associated with the thermal pulses in Sects. 7.3 
to 7.5. 

7.2. Lithium production 

Sackmann & Boothroyd (1992) were the first ones to quan- 
titatively demonstrate how the Cameron & Fowler (1971) 
scenario leads to ^Li synthesis in intermediate-mass AGB 
stars. This major result was obtained by using a time- 
dependent convective diffusion algorithm suited for hot- 
bottom convective envelopes. Indeed, at the base of such 
envelopes, some nuclear reactions involved in ''Li produc- 
tion occur faster than mean turn-over time associated with 
the convective motions (see Fig. ^). 

We also use a time-dependent treatment to study the 
^Li synthesis [see Eqs. (4) in Sect. 2.2.3]. We show in Fig. 
^ the typical resulting ^Be and ^Li profiles inside the con- 
vective envelope of our 6 Mq AGB models with Z = 0.02. 
When HBB occurs (see Sect. 7.1.2), "^He first begins to 
burn through the ^He (^He , 2p) and '^He (a , 7) ^Be reac- 
tions, the second one being significantly slower than the 
first one. As in deep stellar interiors, the ^Be destruction 
by electron capture and proton capture are much longer 
than the time scale associated with the convective mixing, 
the produced ^Be is efficiently up-heaved towards cooler 
regions where it progressively decays in ^Li. Let us em- 
phasize that without fast convective mixing, ^Li could not 
accumulate at the surface. Indeed, at the base of the con- 
vective envelope, (i) ^Be would be partially destroyed by 




r (Ro) 

Fig. 4. Nuclear reaction time scales for the reactions involved 
in the Cameron-Fowler scenario, inside the convective envelope 
of the 6Mq {Z = 0.02) star during the 14th inter-pulse. At 
that time, its base temperature is ~ 75 10® K. r(^Be) and T(^Li) 
correspond to the proton capture on ^Be and ^Li, respectively. 
r("^He) and r("^He-|-Qf) relate to the '^He destruction by ^He and 
*He capture, respectively. Finally, r(^Be-|-e~) is the electron 
capture rate on ^Be. The mean turn-over time of the convective 
motions is also plotted (tconv) 

electron and proton captures and (ii) the ^Li resulting 
from ^Be decay would be quite instantaneously destroyed 
by proton capture too. 

Let us now compare our predictions with Sackmann 
& Boothroyd (1992) results and with observations. A di- 
rect comparison with their Figs. 2a and 2b reveals nearly 
identical features. In Fig. |^, we present the evolution of 
the surface ^Li abundance during the TP- AGB phase of 
all our modeled stars. In our 3 Mq models, the tempera- 
ture at the base of the convective envelope is never high 
enough for ^Be to be produced efficiently, and the sur- 
face ^Li abundance remains almost constant at a very low 
value resulting from the first dredge-up. We also see that 
our 4Mq star with Z = 0.005 begins to produce some 
^Li when it enters the asymptotic regime, but the surface 
abundance of this light element always stays at least one 
order of magnitude below the interstellar medium (or cos- 
mic) value £(^Li) = log ^ + 12 cri 3.3. Note the ^Li destruc- 
tion inside the 4Mq star with Z = 0.02 attesting that its 
convective bottom temperature is high enough to burn ^Li 
through ^Li (p , a) but not enough to produce ^Be by ^He 
burning. Finally, our 5 and 6 Mq models present real high 
^Li surface enhancements, compared to the cosmic value. 
The maximum abundances we obtain range between £(^Li) 
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Fig. 5. ^Be and profiles inside the same convective envelope 
as in Fig. ^ 

~ 4.2 to 5.1, depending on the stellar mass and metallicity. 
Our super-lithium-rich stars are predicted to appear in a 
luminosity range between 2.610* and 3.610*Lq, i.e. for 
Mboi between -6.25 and -6.65. These values are in good 
agreement with those found by Sackmann & Boothroyd 
(1992). 

As the evolution proceeds, the surface ^Li abundances 
finally decrease again very rapidly. This indicates that 
■^He has been almost completely burned in the envelope 
[remember that the '^He ('^He , 2p) is dominant], so that 
the ^Li production can no more be supported and it is 
then destroyed by proton capture. This fact clearly puts 
an upper limit in luminosity to observe super-lithium-rich 
AGB stars. This last feature is also present in Sackmann 
& Boothroyd (1992), even if it appears somewhat later 
during the TP- AGB phase of their most massive objects. 
This could be due to different mass loss rate prescriptions. 

How do these predictions compare with observational 
data (see Fig. ^? In the Magellanic Clouds, the large ma- 
jority of lithium-rich giants are observed within a narrow 
luminosity range, M^oi ~ — 6 to ~ —6.9 (Smith & Lam- 
bert 1989, 1990b; Plez et al. 1993; Smith et al. 1995), 
in very good agreement with the predictions. The ob- 
served abundance range is also well reproduced, even if 
no star has been observed with £(^Li) higher than 4.5 in 
the Clouds. Moreover, the mass estimate from pulsation 
theory for some lithium-rich giants of both Clouds is con- 
sistent with the predicted mass range for surface ^Li en- 
hancement (Smith et al. 1995). Comparisons with ^Li in 
galactic giants are more delicate because of the larger ob- 
servational uncertainties on luminosity and mass. Galactic 



super- lithium- rich AGB stars can exhibit £(^Li) values as 
high as 5.4 (Abia et al. 1991), but these carbon stars have 
Mboi laying between -5 and -6.2, i.e. rather lower than the 
predicted luminosities. This could be due to a bad lumi- 
nosity determination, as it seems difficult to synthesize 
^Li inside the convective envelope of AGB stars with such 
relatively low luminosities. However, most of the galac- 
tic super-lithium-rich AGB stars have lower masses than 
the ones we model here, and they will be addressed in a 
forthcoming paper. 
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Fig. 6. Resulting surface ^Li abundances, in term of e(^Li), 
for the seven computed TP-AGB phases, as a function of the 
stellar surface luminosity. The solid (dotted) lines refer to the 
Z = 0.02 (0.005) models. Observations are shown for bright 
AGB stars of the SMC and LMC (Smith et al. 1995) and of 
the Galaxy (Abia et al. 1991) 



7.3. CNO nucleosynthesis and Fluorine production 

We now review all the species (from C) by indicating in 
which nuclear reactions they are implicated inside the con- 
vective tongue of the successive thermal pulses, as a func- 
tion of the initial mass M and metallicity Z of the com- 
puted TP-AGB models. For the discussion, we refer to 
Figs. 7 to 9 of the electronic version of the paper. The 
mean mass fractions Xa are calculated by taking the av- 
erage of the mass fraction profile of each nuclide a over 
the maximum region covered by the convective tongue of 
each thermal pulse, just after its disappearance, i.e. 

tongue tongue 
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We also refer to Fig. |7| that presents the mean mass 
fractions of key nuclides (calculated as for the convective 
tongue) in each burning region (i.e. the HBS and inter- 
shell regions) compared to their convective envelope con- 
tent, for the 3Mq star with Z = 0.02. Similar figures for 
the other six AGB stars are only available in the elec- 
tronic version. Similarly, Fig. ^ presents the abundance 
variations of the same nuclides inside typical convective 
thermal pulses. Again, similar figures for the other six 
AGB stars can solely be found in the electronic form. All 
these figures are supporting the following nucleosynthesis 
descriptions. 

During the growth of the convective tongue, the mean 
^^C mass fraction first decreases due to the "^He mixing 
from the inter-shell region. Then, at the thermal pulse 
maximum, ^^C is produced by the 3a reaction and also 
enhanced due to the downwards penetration of the convec- 
tive tongue into the He-burnt shell. Although this pene- 
tration is commonly found, its amplitude seems somewhat 
more important in our case, bringing more ^^C into the 
thermal pulses. Note this increased amount of ^^C makes 
the 3DUP episodes more efficient to convert an AGB star 
into a C star. During the first thermal pulses, the amount 
of ^^G that is produced increases due to the higher tem- 
peratures reached at the base of the convective tongue. 
This bottom temperature increase from pulse to pulse is 
a common feature of all the stars and it traduces that the 
thermal pulse intensity increases. The last computed ther- 
mal pulses of the 3Mq object with Z = 0.02 are even hot 
enough to partially destroy ^^C by the ^^C (a , 7) reac- 
tion. 

Most of the engulfed ^^C is destroyed very efficiently 
through -'^^C (a , n) that produces neutrons. The trend of 
the remaining mean ^'^C mass fraction to increase from 
pulse to pulse is due to the fact that stronger thermal 
pulses develop convective tongues reaching regions always 
closer to the HBS, i.e. engulfing more material from the 
inter-shell zone. As mentioned in Sect. 6.1, the bottom of 
the HBS cannot be reached however. 

"'^'^C is only produced by the capture of part of the neu- 
trons on the engulfed ^"^N, through ^"^N (n , p) . Note that 
this reaction contributes to the production of protons. For 
initial masses above 4M0 (3Mq) with Z = 0.02 (0.005) 
however, the convective tongue finally becomes hot enough 
(base temperature > 260 10® K) to significantly destroy 
the produced ^^C by the ^''C {a , 7) reaction. The rela- 
tively low ^''C production inside the 3Mq with Z — 0.02 
object is due to the relatively lower amount of -'^^C that 
is ingested (partially due to the third ^'^C burning shell; 
see Sect. 7.1.1), leading to a smaller quantity of neutrons. 
However, the last thermal pulses of this star ingest increas- 
ing amounts of ^^C and ^"'N from the inter-shell, leading 
to an increasing production of ^'^C. Due to its short life- 
time {ti/2 = 5730 yr), ^"^C partially decays during the 
inter-pulse phases of our less massive AGB models. This 
depletion is however reduced by the "'^^C production inside 



the radiative ^^C burning shell. For the most massive AGB 
stars that have inter-pulse durations shorter that ii/2, 
mostly behaves like a stable nuclide, even if strongly de- 
stroyed by a captures. Consequently, as explained with 
more details in Forestini et al. (1996), intermediate-mass 
AGB stars are globally not important producers of ^''C. 

Apart from the first two thermal pulses (whatever the 
initial mass or metallicity) , is almost completely de- 
stroyed inside thermal pulses through ^^N (a , 7) . This 
reaction, as first noted by Iben (1976), also significantly 
contributes to the energetics of a thermal pulse. Only a 
small part of ^''N is destroyed through the ^"^N (n , p) ^^C 
reaction. Indeed, the inter-shell mean mass fraction of ^''N 
IS rang ed between ~ 10"^ for the 3 Mq star with Z = 0.02 
and 10"^ for the 5M0 star with Z = 0.005, while the 
production level of ^''C is ^ 10^^ (^ 10^®), in mass frac- 
tion, for the Z = 0.02 (0.005) models (except for the 3 Mq 
with Z ~ 0.02 star; see above). 

As are ^^C and ^"^N, ^^N is engulfed in a thermal pulse 
from the inter-shell region. The only significant reaction 
that produces this element inside the convective tongue it- 
self is (p , a) (see below for the ^^0 production). On 
the other hand, ^^N is mainly destroyed by the ^^N (a , 7) 
reaction. Its rate is lower than that of ^^N (a , 7) so that 
-'^^N is efficiently destroyed in all the stars only in the 
asymptotic regime, i.e. when thermal pulses become very 
intense. Remember that it arrives earlier for lower Z ob- 
jects. Note that the 3Mq stars (whatever Z) maintain 
a relatively higher ^^N abundance, especially during the 
full amplitude thermal pulses. Indeed, while these thermal 
pulses efficiently destroy ^^N whatever the initial mass of 
the star, its initial abundance inside the convective tongue 
is higher for 3 Mq objects. This is due to the fact that the 
engulfed inter-shell region of 3Mq stars, that contains 
-'^^N, is significantly thicker (in mass) compared to that of 
more massive ones. 

^^0 being partially destroyed in the HBS, the mixing 
of the inter-shell region reduces its mean mass fraction. 
On the other hand, at the maximum extent of the ther- 
mal pulses, it is enhanced, like ^^C, due to the downwards 
penetration of the convective tongue. From pulse to pulse, 
it also becomes significantly produced by ^^C (0,7) . This 
is especially visible for Z — 0.005 stars, the initial ^®0 con- 
tent of which is lower. This primary ^®0 production by the 
thermal pulses finally stabilizes around ~ 10~^ (in mass 
fraction). The decrease observed during the last thermal 
pulses of the 3M0 AGB star with Z = 0.02 traduces the 
more important role of mixing of inter-shell material (al- 
ready emphasized for ^''N and ^^N). 

The ^^O previously contained in the inter-shell region 
is also mixed in the convective tongue. It is not produced 
inside thermal pulses. On the contrary, it can be destroyed 
mainly through -'^^O {a , n) and ^''O (0,7) , the first re- 
action being in general at least ten times faster than the 
second one. This weakly contributes to produce neutrons 
inside the thermal pulses. This ^^O destruction already op- 
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erates from the first thermal pulses, whatever the initial 
M and Z. fn later thermal pulses, the ^^O mass fraction 
stabilizes around 4 to 810~^°, that represents a balance 
between the amount of ^^O brought from the inter-shell re- 
gion and destroyed by the thermal pulse. Again, our 3 Mq 
models with Z — 0.02 indicate a still dominant role of the 
inter-shell mixing. It is less true for the Z = 0.005 models 
because at lower Z, thermal pulses are somewhat hotter, 
that favors ^^O depletion. 

Contrarily to the case of ^^O, the inter-shell mass frac- 
tion of ^^O is very low 10~^ typically). Its presence 
inside thermal pulse consequently results from its produc- 
tion by the ^"'N {a , 7) reaction. Note that the rapid 
decay to ^^O also significantly contributes to the 
thermal pulse energetics. This element is destroyed by 
three reactions that are, by decreasing order of nuclear 
time scales, ^^O (0,7) , ^^O (a , n) and ^^O (p , a) ^^N . 
This last reaction can even become faster than ^^O {a , n) 
when significant amounts of protons are produced. This 
especially arrives for the late thermal pulses of the 3 Mq 
star with Z = 0.02. The convective tongues being hotter 
for more massive AGB stars and/or later along the TP- 
AGB phase of a given star, the destruction, which is 
very sensitive to temperature, becomes more and more ef- 
ficient from pulse to pulse. This point was already noted 
by Boothroyd & Sackmann (1988b). The thermal pulse 
equilibrium mean mass fraction between production and 
destruction is ranged between 10^^ and 10^^ depending 
on M and Z . 

As for ^*0, ^^F is destroyed in the HBS. It is however 
produced inside thermal pulses. As first demonstrated by 
Forestini et al. (1992), the ^^F synthesis involves many 
nuclear reactions, so that it constitutes a good tracer of 
the nucleosynthesis conditions prevailing inside thermal 
pulses. ^^F is produced by the ^^N [a , 7) reaction. The 
rather rapid increase of the mean ^^F mass fraction dur- 
ing the first thermal pulses of each AGB star is due to 
the increasing importance of the inter-shell ^^N destruc- 
tion. When convective tongues become hot enough to com- 
pletely destroy the engulfed ^^N, the following thermal 
pulses only ingest ^^N from matter that has not expe- 
rienced the preceding thermal pulse (i.e. ^ 50% of the 
convective tongue extension in mass typically; see Ta- 
ble III). However, as we have explained above, ^^N can be 
produced by the chain {a , 7) ^^F (/3+) ^^O (p , a) ^^N . 
Neutrons are mainly coming from the ^'^C (a , n) reac- 
tion, and to a lower extent, from (a , n) , ^^Ne (a , n) , 
^"'^Ne (a , n) (see below) and ^^O (a , n) . Protons are com- 
ing from (n, p) reactions, mostly on ^^Al^^ and ^''N (see 
below). However, (i) only part of the neutrons are avail- 
able to produce protons as many (n, 7) reactions also occur 
(see Sect. 7.5) and (ii) a small (<^ 10%) part of the protons 
are also captured by other nuclides than ^^O [mainly by 
""^^N (p , a) ^^C]. The success of this reaction chain mainly 
depends on the rate of ^^O destruction, mainly by a cap- 
tures (see above) . This destruction becomes more efficient 



from pulse to pulse so that only a small fraction of ^^O is 
available for the (p, a) reaction. On the other hand, when 
convective tongues become very hot (base temperature 
> 28010^ K), i^F begins to be destroyed by ^^F(a,p) 
and ^^F (n , 7) (if significant amounts of neutrons are 
available). In conclusion, three facts concur to consider- 
ably decrease the ^^F production after the first thermal 
pulses (whatever M and Z). Indeed, from pulse to pulse 
and especially for the most massive AGB stars, 

— the amount of ^'^C engulfed in the convective tongue 
of a thermal pulse decreases (due to the partial de- 
struction of ^'^C in the inter-shell region during the 
inter-pulse phase; see Sect. 7.1.1); 

— the amount of ^^N ingested from the inter-shell region 
also decreases; 

— the a captures on ^^O are more frequent; 

— the ^^F destruction becomes significant. This can in 
particular lead to a partial destruction of the ^^F syn- 
thesized in the previous thermal pulses for stars with 
M > 4Mq. 

7.4- NeNa and MgAl nucleosynthesis 

The mean ^"Ne abundance is roughly unchanged in the 
HBS as well as by the thermal pulse nucleosynthesis. In 
number, the contribution of ^^O (a , n) ^"^Ne is rather neg- 
ligible. 

^^Ne is significantly destroyed in the HBS. Inside the 
convective tongues of quite hot thermal pulses (base tem- 
peratures above 27010^ K), it is produced by the small 
proportion of that is destroyed through the ^^O {a , n) 
reaction. There is also a small contribution by "'^^O (a , 7) . 
Only a very small ^^Ne fraction can be destroyed by the 
^^Ne (a , n) at the end of very hot thermal pulses (base 
temperatures above 290 10^ K). 

^^Ne is somewhat produced in the HBS. Its abundance 
is considerably enhanced inside thermal pulses (by typi- 
cally a factor of 100) due to its production by ^^O (0,7) . 
However, its abundance from pulse to pulse becomes con- 
stant (or even slightly decreases inside full amplitudes 
thermal pulses of the most massive stars and/or those 
with lower Z). At their maximum extent, such convec- 
tive tongues are effectively hot enough (base tempera- 
tures > 290 10*^ K) to activate the ^s^e (a , 7) ^'^Mg and 
^^Ne (a , n) ^^Mg reactions. However, as noted below, only 
a small part of the ^^Ne nuclides are converted into ^^Mg 
or ^^Mg, respectively (typically one hundredth). Never- 
theless, this consists in a non-negligible neutron source 
for these thermal pulses, as first noted by Kappeler et al. 
(1990) and emphasized by Straniero et al. (1995). 

Concerning ^■^Na, one has to distinguish between the 3, 
4, 5 M0 with Z = 0.02 and 3 Mq with Z = 0.005 stars on 
the one hand, and the 6Mq with Z = 0.02, 4 and 5M0 
with Z = 0.005 stars on the other. Inside the former ones, 
^^Na is significantly produced at the bottom of the HBS 
and very slightly produced by the HeBS [during the inter- 
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pulse phase, through the ^°Ne (a , p) reaction] . Further- 
more, whatever M and Z, ^^Na is almost unchanged by 
the thermal pulse nucleosynthesis. As a consequence, its 
mean mass fraction inside the convective tongue does not 
evolve significantly from pulse to pulse. At the opposite, 
inside the latter ones, the bottom of the HBS is hot enough 
to somewhat destroy ^^Na while it is produced in its up- 
per part. However, the HeBS substantially produces this 
nuclide during the inter-pulse phase. Consequently, when 
the convective tongues of such stars penetrate the inter- 
shell region, they engulf material that is impoverished in 
^^Na; this explains its mean mass fraction decrease from 
pulse to pulse. 

The ^"'Mg abundance is almost unchanged by the ther- 
mal pulse nucleosynthesis too. Furthermore, the HBS of 
SMq stars is not hot enough to destroy it significantly. 
Its mean mass fraction in the convective tongue region 
of such stars consequently remains almost constant from 
pulse to pulse. The situation is quite different for more 
massive stars (whatever Z) as they destroy ^^Mg at the 
bottom of their HBS. Dilution of inter-shell material by 
the successive convective tongues thus traduces in a mean 
mass fraction decrease with time. 

Intermediate nuclide of the MgAl chain, ^^Mg is mod- 
ified by the proton burning inside the HBS. We refer to 
Sect. 7.1.1 for more details (dependence on M and Z). 
However, the evolution of its mean mass fraction inside 
the convective tongues of thermal pulses is dominated by 
its production through the ^^Ne (a , n) reaction [operat- 
ing at a quite comparable rate than ^^Ne (a , 7) , following 
our present knowledge of the corresponding nuclear reac- 
tion rates]. Indeed, even if this reaction is very slow inside 
thermal pulses of relatively low mass AGB stars (espe- 
cially for their first thermal pulses), ^^Ne is much more 
abundant than ^■''Mg (by typically a factor of one hun- 
dred) in that region, that explains the sensible increase of 
the ^^Mg mean mass fraction. 

^^Mg is slightly (significantly) depleted in the HBS 
of the Z = 0.02 (0.005) models. On the other hand, 
it is substantially produced by the thermal pulse nucle- 
osynthesis (i) by the ^^Ne(Q!,7) reaction (same remark 
as for ^^Mg above) and, to a lower extent, (ii) by the 
^^Al^ (n , p) reaction (see below). The slope with which its 
mean mass fraction increases form pulse to pulse is conse- 
quently higher for hotter convective tongues, i.e. steeper 
for full amplitude thermal pulses of massive AGB stars. 

^^Al^ is one of the most important by-product of the 
HBS (see above). With a lifetime fi/2 = 7 10'' yr, i.e. much 
longer than the inter-pulse duration of intcrmcdiatc-mass 
AGB stars, it essentially behaves like a stable nuclide in 
the inter-shell region. When engulfed inside the convective 
tongue of a thermal pulse, it is mainly destroyed through 
the ^^Ais (n , p) reaction and, at least ten times slower, 
by ^^AF (n , a) . As a result, ^^Al^ is the most important 
proton source inside thermal pulses, whatever M and Z. 
Its destruction is however partial, due to the relatively 



low neutron abundance (see Sect. 7.5 below). Note that 
^^AF (n , p) is also the principal neutron source in the 
radiative ^^C burning shell, as reported by Wasserburg et 
al. (1994) too. 

^^Al is significantly produced in the HBS, especially 
in the most massive AGB stars. It docs not significantly 
participate to the thermal pulse nucleosynthesis. As a con- 
sequence, it is engulfed by the convective tongue from the 
inter-shell region and its corresponding mean mass frac- 
tion slightly increases from pulse to pulse. Note however 
that with time, most of the H-burning of the 5 Mq AGB 
star with Z = 0.005 occurs inside its convective envelope. 
This considerably reduces the ^''Al production. 

Si, P and S are almost not concerned by charged parti- 
cle nuclear reactions inside the convective tongue of ther- 
mal pulses. On the other hand, they are also rather un- 
changed in the HBS. The only role played by these ele- 
ments concerns neutron captures that we now briefly dis- 
cuss. 

7.5. Neutron production 

Many evolved AGB stars show considerable overabun- 
dances in elements liciavier than iron (sec e.g. Smith & 
Lambert 1986). The original discovery by Merrill (1952) of 
Tc in the spectra of some S stars definitively proved that 
these heavy elements are synthesized inside AGB stars. 
Iben (1975) first demonstrated how such species can in- 
deed be produced inside the thermal pulses by slow neu- 
tron captures (the so-called s— process), those neutrons 
coming from the ^^Ne (a , n) reaction. However, Iben's 
computations concerned a very massive AGB star (7 Mq), 
while ever since, all the observed stars showing s elements 
enrichment have been identified as being of lower mass 
(i.e. < 6 Mq and for many of them < 3 Mq). More recent 
observations (Aaronson & Mould 1985) and evolution- 
ary models (Malaney & Boothroyd 1987) then strongly 
suggested that the major neutron source had to be the 
"'^^C (a , n) reaction. 

^^C is engulfed by the convective tongues of thermal 
pulses from the inter-shell region where the HBS ashes ac- 
cumulate. However, the amount of ^"^C spread out by the 
HBS is by far (i.e. at least a factor of 10 or more) insufii- 
cient to explain the s— process that requires more than one 
neutron by seed ^^Fe. This so-called "s— process mystery" 
(Sackmann & Boothroyd 1991) is a common failure of all 
the recent stellar evolution models of TP- AGB stars. It is 
moreover reinforced by the fact that other lighter nuclides 
also capture part of the available neutrons (mainly, by in- 
creasing order of atomic mass, ^^C, ^^N, ^''O, ^'^F, ^%Ig, 
25Mg, 26Alg, 28si^ 29si^ sogi^ 3ip^ 323^ ssg^ _ _ _ )_ Que has 

to note that while the neutron abundance has to be in- 
creased to solve the s— process problem, the present com- 
putations probably slightly overestimate the abundances 
of the above mentioned nuclides. 
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More specifically, in our thermal pulse models, the ra- 
tio of the engulfed mass of ^^C to the mass of ^^Fe in- 
side the convective tongue at its maximum extent ranges 
between 0.002 to 0.03. This upper hmit is reached dur- 
ing the asymptotic regime of all our modeled stars and 
thus appears to be rather independent of the initial mass 
or metallicity. Inside their first thermal pulses however, 
3 Mq stars already produce neutrons (ratio roughly equal 
to 0.02) while at the opposite, the corresponding thermal 
pulses of 6 Mq stars do not engulf substantial amounts of 
^^C. Furthermore, we noted (see Sect. 7.4) that the ^^Ne 
neutron source operates after the maximum of very in- 
tense thermal pulses, at a rate that is growing with the 
maximum temperature of the convective tongue, i.e. with 
initial mass and/or pulse number. 

Lot us stress that such low quantities of ^''C ingested 
by thermal pulses are partly explained by the partial ^^C 
radiative burning, operating at the bottom of the inter- 
shell region during the inter-pulse phase (see Sect. 7.1.1). 
Consequently, it seems clear that the s— process occuring 
inside TP-AGB stars has two distinct origins (in time and 
space): 

— during the inter-pulse phase, radiative s— process oc- 
curs in a radiative zone that is the lower part of the 
inter-shell region; the only neutron source is ^^C; 

— during and inside the convective tongue associated 
with thermal pulses, the classical s— process takes 
place with two possible neutron sources: (i) the re- 
mainder amount of ^^C that is engulfed by the con- 
vective tongue and, to a lower extent (maximum a few 
percents), (ii) the ^^Ne, close to the end of hot thermal 
pulses. Let us emphasize that the (a, n) and accom- 
panying (n, 7) reactions have non-negligible energetic 
effects on the structural evolution of the convective 
tongues. 

We compute again some thermal pulses by treating 
together nucleosynthesis and time-dependent convective 
mixing through Eqs. (4). We observed an abundance gra- 
dient inside the convective tongue for a few nuclides that 
are involved in very fast nuclear reactions (neutrons, pro- 
tons and some unstable nuclides like ^*F or ^^Na). Glob- 
ally however, the final abundances resulting from the ther- 
mal pulse nucleosynthesis were similar to within ^ 10% 
compared to standard computations. More specifically, 
neutrons were systematically found to be much more 
abundant (by typically five orders of magnitude) at the 
base of the convective tongues, meaning that the ingested 
^•^C is first transported down before it is destroyed by 
^•^C (a , n) . The liberated neutrons are very rapidly cap- 
tured, as attested e.g. by a greater proton abundance at 
the base of the convective tongue too. This leads us to 
the conclusion that matter that is actually irradiated by 
the neutron flux is only the bottom part of each convective 
tongue. The number of free neutrons by seed ^^Fe is so en- 
hanced. Such an approach has already been suggested by 



Malaney et al. (1988). The conclusion, however, remains 
that in order to reproduce the observed distribution of 
s elements in the primitive solar system, higher neutron 
fluxes are still needed. 

The most natural way to conciliate observations and 
theoretical models concerning the production of the s el- 
ements should be to increase the amount of '^"^C in the 
inter-shell region. This can be done if protons are trans- 
ported down, e.g. during the 3DUP. ^^C could then be pro- 
duced by the ^"^C (p , 7) reaction. Such a scenario would 
enhance the s— process in both sites where it can occur. In 
our opinion, the present failure to build this extra-amount 
of ^'^C could be related to our bad treatment of convective 
boundaries, especially when they rapidly penetrates very 
inhomogeneous regions, like it is the case during the 3DUP 
events. The possible occurrence of such a slow-particle 
transport process able to transport protons down into the 
inter-shell region and its exact efficiency will be quantita- 
tively investigated in the context of low-mass AGB stars. 
These stars have indeed longer inter-pulse and 3DUP du- 
rations that could allow to transport signiflcant amounts 
of hydrogen downwards (see Sect. 9). 

Straniero et al. (1995) have recently shown that such 
enhanced amounts of ^^C in the inter-shell region could 
indeed allow the radiative s— process to occur. Most in- 
terestingly, they found that the resulting signatures are 
rather similar to those coming from classical s— process 
computations and are even reached faster. This work is 
very important as it clearly demonstrates that this so- 
called "^^C pocket" allows to explain the observed s ele- 
ment synthesis in AGB stars. 

7. 6. Resulting abundances of each nuclide inside the vari- 
ous nucleosynthesis sites 

In order to summarize our detailed discussion of nucle- 
osynthesis in the previous sections, we now identify wliic;li 
specific nuclear region, inside a TP-AGB star, mainly con- 
tributes to the surface abundance change of each nuclide, 
when the third dredge-up occurs. 

Concerning the light elements, let us just recall (see 
Sect. 7.2) that ''Li can be produced by HBB in the con- 
vective envelope of > 4.5 Mq (4Mq) AGB stars with 
Z = 0.02 (0.005). 

When a third dredge-up occurs, the convective enve- 
lope first penetrates the inter-shell region that is depleted 
in ^^C and then reaches a deeper region previously en- 
riched in ^^C by the convective tongue. As a result, ^^C 
globally increases in the convective envelope. The situa- 
tion is reversed for ^'^C. After a 3DUP, the envelope is 
somewhat impoverished in ^^C as this nuclide is almost 
completely destroyed inside the convective tongue. Conse- 
quently, the isotopic ratio ^^C/^^C significantly increases 
at each 3DUP, all the more as the convective envelope 
mass is reduced. However, in our > 5Mq (4Mq) TP- 
AGB stars with Z = 0.02 (0.005), HBB is strong enough 
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to partially convert ^^C in ^^C through the CN cycle so 
that the surface ^^C/^'^C drastically decreases. For these 
objects, the 3DUP and HBB have opposite effects. Last 
but not least, -^^C, somewhat produced by thermal pulses 
and the third burning shell, is up-heaved to the surface 
during 3DUP events, where it however substantially de- 
cays in the less massive AGB stars. 

"^^N is produced in the inter-shell region but is strongly 
destroyed inside thermal pulses. ^^N is destroyed in the 
intcr-sliell region while it is somewhat produced by the 
first few thermal pulses and then destroyed inside full 
amplitude thermal pulses, i.e. those that can be followed 
by a 3DUP. As a consequence, the surface isotopic ratio 
^''N/^^N is expected to slightly increase from 3DUP to 
3DUP. Again, as the CN cycle operates inside the con- 
vective envelope of the most massive AGB stars, ^^N is 
produced from ^^C and ^^N is destroyed, thus leading to 
very high ^''N/^^N ratios. If HBB occurs at very high tem- 
peratures (i.e. ^ 10^ K) a rather long time enough (i.e. 
for the lowest mass loss rate), the approach towards the 
CNO bi-cycle equilibrium leads to an increase of the "^^N 
abundance and consequently, a decreasing ^"'N/^^N ratio. 
This mainly concerns our 6M0 (5Mq) AGB star with 
Z = 0.02 (0.005). 

-'^^O, like ^^C, is depleted in the inter-shell region, but is 
produced inside the thermal pulses. During 3DUP events, 
its surface abundance very slightly increases. ^'^O is always 
considerably depleted in both the inter-shell and thermal 
pulse regions, so that a 3DUP decreases its surface abun- 
dance. Finally, at the end of full amplitude thermal pulses, 
^^O has a lower abundance as inside the convective enve- 
lope. As it is also completely destroyed in the inter-shell 
region, its abundance drastically decreases from 3DUP to 
3DUP. Consequently, both ^sq/^o and i^o/isq surface 
isotopic ratios significantly increase during each dredge- 
up. If HBB is strong enough however, the operation of 
the ON cycle increases the amount of ^^O to the ex- 
pense of and substantially destroys ^^O. So, in massive 
AGB stars, i^o/i^Q can decrease while ^^Of^^O further 
increases much more than in stars that do not experience 
HBB. 

is destroyed in the HBS but it is significantly pro- 
duced by thermal pulses. As a consequence, 3DUP events 
must lead to correlated ^^F and ^^C surface enhancements, 
since both elements are produced in the same region. 

As already explained, ^°Ne is almost unchanged in 
the HBS and thermal pulse nucleosynthesis. ^^Ne is de- 
stroyed inside the inter-shell region but produced inside 
thermal pulses, this production slightly increasing from 
pulse to pulse, whatever M or Z. Third dredges- up con- 
sequently slightly decreases the ^°Ne/^^Ne isotopic ratio. 
In our 6M0 (SMq) AGB models with Z = 0.02 (0.005) 
however, ^^Ne is partially destroyed by HBB. On the other 
hand, the ^°Ne/^^Ne ratio substantially decreases from 
3DUP to 3DUP due to large ^^^Ne production in the HBS 
and especially inside the thermal pulses. 



During 3DUP events, the surface abundance of ^^Na 
slightly increases for AGB stars with initial masses up to 
4M0 (3 M0) and Z = 0.02 (0.005), due to its production 
in the HBS. The hotter HBS of more massive AGB stars 
do not produce ^^Na, globally. Our 6Mq (5Mq) models 
for the same respective Z even slightly destroy ^'^Na by 
HBB. 

^"'Mg is somewhat destroyed in the HBS of our most 
massive AGB stars, especially those with Z = 0.005, while 
its abundance does not change significantly after the ther- 
mal pulse nucleosynthesis. ^^Mg is depleted in the inter- 
shell region of our 3 Mq models with Z = 0.02 as well as 
in the Z = 0.005 models of all masses. It is however pro- 
duced in the HBS of > 4M0 AGB stars with Z = 0.02. 
Nevertheless, its surface abundance evolution following 
3DUP events is mainly conditioned by its large abun- 
dance (compared to the envelope one) in the region where 
thermal pulses occur, especially inside the most massive 
AGB stars. Therefore, the surface ^^Mg/^^Mg isotopic ra- 
tio necessarily decreases in all AGB stars. The same is 
true for the ^''Mg/^^Mg ratio that roughly behaves in a 
similar way. Reasons are also rather similar, even if ^^Mg 
is less modified than ^^Mg in the HBS. At the surface of 
our 6Mq (5M0) with Z = 0.02 (0.005) AGB stars, the 
abundances of both ^^Mg and ^^Mg can be enhanced to 
the expense of ^^Mg, if HBB is strong enough, leading to 
very low ^%Ig/^^Mg and ^^Mg/^®Mg ratios. As we will see 
in Sect. 8.1 however, this only concerns the end of the TP- 
AGB phase in extreme situations allowing HBB to occur 
a long time enough (e.g. in case of relatively low mass loss 
rates). 

Only a very small amount of ^^Al^ is present in the 
convective envelope of E-AGB stars as a result of the first 
(and eventually the second) dredge(s)-up. Along the TP- 
AGB phase, ^^Al^ is largely produced in the HBS and most 
of it accumulates in the inter-shell region. Inside full am- 
plitude thermal pulses however, it is somewhat destroyed. 
For its part, ^^Al is slightly produced in the HBS, espe- 
cially inside the most massive AGB stars. Consequently, 
^^AF/^^Al significantly increases with time during the 
TP-AGB phase, due to the repetitive 3DUP episodes. Of 
course, dilution by convective mixing being more impor- 
tant for increasing M, the highest ratios are expected at 
the surface of the less massive AGB stars. However, the 
same extreme cases mentioned in the case of the Mg iso- 
topic ratios can lead to very high ^^Al^/^'^Al ratios for the 
most massive AGB stars. 

^From Si, the 3DUP does not change significantly the 
surface isotopic ratios. However, as stressed in Sect. 7.5, 
when models will be able to reproduce the s element en- 
hancements at the surface of evolved AGB stars. Si and 
S isotopic ratios could be slightly modified by the high 
neutron flux. 
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Fig. 8. Variation of the mean mass fractions of the same 20 nuchdes as in Fig. |7|, due to the thermal pulse nucleosynthesis, for 
the 3Mq {Z = 0.02) star, inside an early thermal pulse (the 3d) and an asymptotic one (the 18th). Arrows indicate the mean 
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8. Wind composition along TP-AGB phase 

In this last section, we present our predictions concerning 

(i) the evolution of all the surface isotopic ratios along the 
TP-AGB phase (by including our evolutionary models and 
extrapolations up to the convective envelope removal) and 

(ii) the total amount of material ejected at different ages 
for all the species up to ^^Si. We refer to the previous 
sections for detailed explanations. Thus, in the following, 
we will just comment these predictions. 

8.1. Surface isotopic ratios along TP-AGB 

Figs. I and |l^ display, for the SMq AGB star with 
Z = 0.02, the evolution of the isotopic ratios as a function 
of the total remaining mass that decreases due to mass 
loss. Similar figures for the other six AGB stars are only 
displayed in the electronic version. 

Let us stress that we stop our calculations at the last 
thermal pulse that occurs before the convective envelope 
mass is reduced below 0.2 Mq. This stage indeed roughly 
corresponds to the PN ejection that we do not maintain to 
model here. This explains the quite large remaining total 
stellar masses in the figures, especially for those computed 
with the highest mass loss rates. 

The first comment concerns the formation of carbon 
stars. 3M0 AGB stars with Z = 0.02 and 3 and 4Mq 
AGB stars with Z = 0.005 become carbon stars dur- 



ing their TP-AGB phase. The surface becomes 
greater than unity more rapidly for (i) decreasing total 
mass and/or (ii) decreasing Z. This metallicity depen- 
dence is confirmed by the increase of C stars compared to 
M stars from the galactic bulge to LMC and SMC (Blanco 
et al. 1978). Our 6Mq (SMq) AGB stars with Z = 0.02 
(0.005) can also become C stars with our lowest mass loss 
rates. In such cases indeed, in spite of very efficient HBB 
transforming ^^C in ^'^C and ^''N, many 3DUP events can 
still occur with very reduced envelope masses, so that the 
relatively short inter-pulse phases are not long enough to 
substantially destroy ^^C. We will see however, together 
with other isotopic ratios, that these very low mass loss 
rates seem actually extreme. Let us emphasize that the 
only super-lithium-rich star (see Sect. 7.2) we found to be 
also a carbon star is the 5Mq with Z = 0.005 one (with 
our medium mass loss rate). However, as already men- 
tioned in Sect. 6.2, observations indicate that Z — 0.02 
carbon stars are formed easier (i.e. at lower luminosities) 
than predicted by all the models. Let us also stress the 
very high ^^C/^'^C isotopic ratios that can be expected 
at the surface of < 4Mq AGB stars, whatever Z. On 
the other hand, more massive AGB stars experience HBB 
that progressively leads the ^^C/^^C isotopic ratio close 
to its CN cycle equilibrium value (formation of J stars). 
The case of 4Mq stars with Z = 0.005 very clearly shows 
the competition between the 3DUP episodes (increasing 
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Fig. 7. Mean mass fractions of the 20 most relevant nuclides 
inside the convective envelope, the HBS and the inter-shell 
region, for the 3Mq {Z = 0.02) star 



i^C/i'^C) and HBB (decreasing it). For lower mass loss 
rates, the 3DUP events are more numerous for the same 
total mass decrease, leading to a more rapid i^C/^'^C in- 
crease. However, the corresponding longer TP- AGB phase 
allows the base of the convective envelope to reach higher 
temperatures that finally destroy ^^C through the CN cy- 
cle. This efficient transformation of ^^C to ^^C also reduces 
the i^^C/i^O ratio, i.e. the possibility to form carbon stars. 
It happens if the temperatures at the base of the con- 
vective envelopes become higher than ~ 75 10^ K. Such 
stars have corresponding surface luminosities > 25000 Lq 
(30000 Lq) for Z = 0.02 (0.005), i.e. Ahoi < -6.25 (- 
6.45). We consequently predict that C stars cannot be 
formed at higher luminosities than this threshold, due to 
strong HBB. This is in agreement with Boothroyd et al. 
(1993) predictions. This also corresponds to the observed 
threshold for AGB stars of the LMC (Cohen et al. 1981). 

Finally, note the surface pollution in ^^G that is more 
important, compared to that of ^^C, for decreasing total 
masses. For 3Mq AGB stars, the surface ^^C/^'^C ratio 
is quite independent of the mass loss rate. This is due to 
the inter-pulse duration that is significantly greater than 
the ^'^C lifetime, so that the ^'^C decay is the dominant 



factor for the ^^C/^^G ratio evolution at the surface of 
such stars. The quite complicated behavior of the ^^C/"'^''C 
isotopic ratio at the surface of our most massive models 
is again due to the competition between the 3DUP that 
adds ^^C into the convective envelope and HBB that can 
partially destroy it. By extrapolation, it is reasonable to 
expect even lower ^^C/^^C ratios in low-mass AGB stars. 
As mentioned in Forestini et al. (1996), "'^^C could become 
detectable for < 3 Mq evolved AGB stars. 

Due to HBB, those stars showing very low surface 
12(--;^13q ratios are also expected to have very high ^'^'N/^^N 
ones. It has to be stressed indeed that in our most massive 
AGB stars, the CN cycle operating inside the convective 
envelope produces large amounts of ^'^N. 

The I'^O/i'^O and ^'^O/^^O surface isotopic ratios also 
evolve as expected from our discussions in Sect. 7. Note 
that quite low ^^O/^'^Q and very high ^^O/^^O ratios, 
which are clear signatures of the operation of the ON 
cycle at the base of the convective envelope (HBB), ap- 
pear rather simultaneously with the CN cycle signatures 
in our most massive objects. Also remark the progres- 
sive I'^O/i^O decrease in our 6 Mq (5 Mq) AGB star with 
Z = 0.02 (0.005). For these objects, the temperature at 
the base of the convective envelope finally becomes high 
enough for ^^O to partially burn. 

The ^^F surface enhancements we predict are displayed 
in a special form, corresponding to the observation pre- 
sentations adopted by Jorissen et al. (1992). Compared to 
their Fig. 8, two conclusions can be stated. 

— There is a clear observed correlation between the ^^F 
excess and the degree of ^^C surface pollution, that 
is perfectly understandable from the identified nucle- 
osynthetic processes and the mixing episodes (see Sect. 
7). Furthermore, our models clearly indicate a trend: 
the less massive the AGB star, the more efficient the 
i^F surface enhancement. In fact, our most massive 
models even show a surface ^^F depletion due to its 
partial destruction by HBB. 

— Our Z = 0.02 models well reproduce only the less 
i^F enriched stars in the versus i^C/isq di- 
agram. Based on the trend just mentioned, we predict 
that the AGB stars showing very large ^^F surface en- 
hancements have initial total masses < 3 Mq . We plan 
to verify this statement, i.e. that low-mass AGB stars 
are the most efficient ^^F galactic producers (see Sect. 
9). 

As expected, the ^^Mg/^^Mg and 24Mg/26Mg sur- 
face isotopic ratios are much less sensitive to the 3DUP 
episodes and HBB. More precisely, only our most massive 
models with the lowest mass loss rates show large ^''Mg 
depletions due to a very strong HBB. Again, these are 
however quite extreme situations. 

Last but not least, we predict substantial ^^Al^ sur- 
face enhancements due to the repetitive dredges-up of the 
HBS and inter-shell regions. The 26A1S/27A1 isotopic ratio 
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tends to decrease with increasing total mass due to (i) the 
greater dilution as mixing occurs in a more massive con- 
vective envelope and (ii) the thinner HBS and inter-shell 
regions (in mass). Nevertheless, ^^Al^/^'^Al could well be 
increasing again if ^^AF was produced inside the convec- 
tive envelope itself by strong HBB. Note however that such 
a very high ^^Al^ production level is actually expected to 
be accompanied by very low Mg isotopic ratios, as ^^Al^ 
is produced by the MgAl chain in the extreme situations 
mentioned above. Let us finally mention that Wasserburg 
et al. (1994) presented computations in which the evo- 
lution of the ^^AF surface abundance has been followed 
self-consistently, by including the effect of neutron irradi- 
ation due to an enhanced amount of ^^C in the inter-shell 
region (to engender a s— process). 

Our results compare very well with those obtained by 
Boothroyd et al. (1994, 1995) concerning the O isotopic 
ratios in intermediate-mass stars. They also mention the 
problem to reproduce the ^^O/^^O ratio, especially in low- 
mass carbon AGB stars. We cannot go further in this com- 
parison without having yet computed the TP- AGB phase 
of low-mass stars. Boothroyd et al. (1995) already sug- 
gested that to conciliate all the data with observations of 
low-mass stars, one probably has to take into account of 
slow-particle transports inside radiative zones, what they 
called "cool-bottom precessing" . That is what we plan to 
do (see Sect. 9). 

8.2. Other comparisons with observations 

Evolutionary models of TP-AGB stars can be strongly 
constrained by the confrontation between predicted and 
observed surface isotopic ratios, especially if different ra- 
tios are observed for the same stars. For most of the ob- 
served AGB stars, only the C and O isotopic ratios have 
been determined (see e.g. Harris & Lambert 1984; Harris 
et al. 1985; Lambert et al. 1986; Harris et al. 1987; Harris 
& Lambert 1987; Kahane et al. 1992). Globally, we found 
agreements. General features come out. 

— As the temperatures required at the base of the con- 
vective envelope to allow ''Li production (see Sect. 7.2) 
are significantly lower than those allowing efficient ^^C 
destruction by the CN cycle, we predict that lumi- 
nous C stars enriched in ^Li can exist, especially for 
lower Z objects (that have lower ^^O content). There 
are indeed a few % of the observed C stars that are 
also super-lithium-rich and there are more numerous 
in the LMC compared to our galaxy (see e.g. Smith 
& Lambert 1990b for the LMC and Feast 1974 for the 
Galaxy). 

— There is no observed C star with ^'^C/^^C < 15. Most 
of the J stars show ^^C/^^C ratios between 3 et 4, i.e. 
very close to the CN cycle equilibrium value. 

— J stars also show rather low ^^0/^''0 isotopic ra- 
tios while only upper limits have been found for the 
16Q/18Q jyiQg^ q£ ^jjg Q^i^gj- g^g^j-g present 



16q^i7q ratios ranged between 300 and 4000 and 
16Q/180 ratios ranged between 500 and 5000. Only 
stars with rather high ^^0/^''0 ratios are difficult to 
explain with our intermediate-mass AGB models (see 
also Boothroyd et al. 1995). Also note that the cor- 
responding O isotopic ratios in the ISM are ^^O/^^O 
= 2460;^^^ and ^^O/^^O = 675t2oo- Again, there is a 
problem with ^^O. Let us however stress that impor- 
tant uncertainties remain concerning the nuclear cross 
sections of proton captures on ^''O. On the other hand, 
low-mass AGB stars are present among those observa- 
tions. 

Presently, there is only one star for which all the iso- 
topic ratios have been determined. It is Cw Leo. This C 
star is one of the nearest evolved AGB star ever observed. 
Its distance is estimated to d — 150 ± 50 pc (see e.g. 
Claussen et al. 1987; Keady et al. 1988; Griffin 1990). This 
corresponds to a surface luminosity L ~ 20000 ± 5000 Lq 
and a present mass loss rate of roughly (2.5 ± 1.0) lO^"' 
Mq yr"""^. Such a high wind makes it quite optically thick. 
A very extended circumstellar envelope has been detected 
to surround Cw Leo, IRC+10216. IRC+10216 shows evi- 
dences for a bi-polar structure (see e.g. GueHn et al. 1993), 
attesting that Cw Leo is very evolved and close to the end 
of its TP-AGB phase. Consequently, its present convec- 
tive envelope mass is probably very reduced 0.5 Mq 
typically) . 

The circumstellar envelope of Cw Leo allowed very ac- 
curate determinations of the isotopic ratios (see Kahane 
et al. 1988 and 1992; Guehn et al. 1995; Forestini et al. 
1996): 12C/13C = 44+1, ^^Cf^^C > 1400, ^^N/^^N > 5300, 
160/170 = MOtlfo, ^^0^0 = 1260tl\l, ^^Mg/25Mg 
- 7.6^?:?, 24Mg/26Mg = 6.5t^;^ and ^'^AV/^^Al ^ 0.01 
(to be confirmed). Furthermore, the Si isotopic ratios ap- 
pear to have interstellar medium values, i.e. unaltered by 
the star. As derived in Guclin et al. (1995) and Forestini et 
al. (1996), almost all these isotopic ratios seem to be con- 
sistent (within the error bars), following our predictions, 
with an AGB star that had a main sequence total mass 
between 4 and 5Mq. However, i^O/i^O is again not well 
explained. As various surface isotopic ratios (especially 
i^N/i^N and i^O/i^O) are significantly changing between 
4 and 5 Mq AGB models, we are presently computing the 
TP-AGB phase of a 4.5 Mq star with Z = 0.02 in order 
to improve this last prediction. 

Finally, another very interesting way to constraint stel- 
lar evolution models of AGB stars is to compare the pre- 
dicted surface isotopic ratios with those determined in 
primitive meteorites (mainly through grains included in- 
side carbon rich chondrites). Most of the i^C/^C, I'^N/i^N 
and ^^Al^/^^Al isotopic anomalies measured in grains (see 
e.g. Ott 1991; Zinner et al. 1991a) are consistent with 
those coming from the surface of evolved TP-AGB stars 
of intermediate- mass. This is of course in favor of the pres- 
ence, when the solar system has been formed, of one or a 
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few AGB stars in its surroundings. Such grains are indeed 
well known to be formed in the cool atmosphere of evolved 
AGB stars (i.e. when mass loss rates become rather high). 
Concerning the measured O isotopic ratios, the situation 
is not so clear as many grains have ratios significantly 
out of the corresponding predicted ranges. Some of these 
grains at least can have been formed in low-mass AGB 
stars. Again, as mentioned in Sect. 8.1, Boothroyd et al. 
(1995) recently suggested that in such stars, cool-bottom 
processing can help to reduce these discrepancies. On the 
contrary, some other grains perfectly agree with our pre- 



dictions. This is for example the case of a ~ 3 /im AI2O3 
grain of the Bishunpur LL3.1 chondrite by Huss et al. 
(1994). The derived ratios, i^o/i^o = 385 ± 5, ^^O/^^O 
= 853±30 and ^^AF/^^Al = (1.7±0.2) 10"^, are in perfect 
agreement with ratios corresponding to a 3— 4 Mq AGB 
star of nearly solar metallicity. Finally, other SiC ^^Al^- 
rich grains raise problems, like that discovered by Zinner 
et al. (1991b), exhibiting a very high ^'^AF/^-'Al ratio of 
~ 0.23. It seems hard to explain such a value from our 
models. However, this grain also exhibits an anomalously 
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large *'*Ca/^°Ca isotopic ratio that seems to exclude its 
AGB star origin. 

8.3. Yields 

Tables | to |l^ give the total mass ejected by the wind 
(in Mq) for various elements (i) aX the top of the first 
RGB, (a) at the beginning of the TP- AGB phase and 
(iii) at the end of the TP- AGB phase (thus including our 
extrapolated nucleosynthesis calculations). As such data 
mainly concerns the chemical evolution of galaxies, we 
also present the net yield (positive or negative) for each 
element. The data of Tables H to 10 correspond to the 



medium mass loss rates (see Sect. 8.1), i.e. the solid lines 
of Figs. H and |l^. Brief comments are made about these 
predictions in the electronic version. 

9. Conclusions and prospects 

We tried to confirm and improve our present knowl- 
edge of the structural and nucleosynthetic evolution of 
intermediate-mass AGB stars. To do that, we have first 
presented in detail (i) the physics and numerical aspects 
of our stellar evolution code, (ii) results concerning evo- 
lution phases prior to the AGB one and (in) comparisons 
with other works. These informations are essential to bet- 
ter understand and appraise the large set of predictions we 
discuss about the thermally pulsing AGB stage. We also 
compare our predictions with various observations (at dif- 
ferent evolutionary stages) in order to clearly identify the 
processes that should be included in future models. 

Let us just mention global key features stemming from 
our intermediate-mass evolution models. 

— Concerning the structure (intensities, temperatures) of 
the thermal pulses as well as the thermal properties 
at the base of the convective envelope, the most deci- 
sive quantities are the total mass and core mass. Both 
change with time in a way that is governed by the in- 
creasing mass loss rate. This translates into a great 
sensitivity of the various nucleosynthesis processes to 
the adopted mass loss rate. Consequently, the evolu- 
tion of the thermal pulse characteristics crucially de- 
pend on this last quantity, unfortunately still rather 
badly known. In addition, the lower the mass loss 
rates, the longer the TP-AGB phase, the more third 
dredge-up events have time to occur and the greater 
are the surface abundance changes due to HBB. This is 
clearly demonstrated by the high sensitiveness of the 
predicted surface evolution of almost all the isotopic 
ratios to the mass loss rate dependence on time. 

— On the other hand, at a given total stellar mass, core 
masses are greater for lower metallicities, as a result 
of the central burning phases. As a consequence, for 
the reason mentioned at the former point, TP-AGB 
stars with Z = 0.005 mostly behave like Z = 0.02 ones 
with a total mass greater by ~ 1 Mq. However, due to 



important differences in initial compositions that also 
influence the nuclear reactions, this is not completely 
true as far as isotopic ratios are concerned. At the end 
of the TP-AGB phase indeed, some yields are quite 
different for both Z we have studied. 

— The neutron captures on intermediate-mass and heavy 
nuclides occur in two distinct sites inside TP-AGB 
stars, namely (i) at the base of the inter-shell zone dur- 
ing the inter-pulse phase and, as already known, (ii) 
inside the convective tongues associated with thermal 
pulses. Along the asymptotic TP-AGB, the amount of 
produced neutrons is rather independent of the initial 
total mass and metallicity (at least for intermediate- 
mass objects). 

— Concerning chemi- 
cal evolution of galaxies, intermediate-mass stars, at 
the end of their existence, appear to be rather signifi- 
cant producers of ^He, ^He, ^Li (most massive and/or 
higher Z stars), "C, "N, ^^O, ^^F, ^^Ne, ^^Na, ^^Mg 
and ^^Al. A found trend, compared to observations, 
however clearly supports the idea that most of the ^^F 
present in the interstellar medium has to come from 
low-mass AGB stars (that point has to be confirmed 
by forthcoming evolutionary models). 
Intermediate-mass AGB stars also substantially pro- 
duce radio- nuclides, namely ^^C and ^^Al. The former 
one mostly comes from lower mass stars while it is the 
contrary for the last one. 

These stars partially deplete the interstellar medium 
content in ''Li (less massive stars), ^^N, ^^O, ^*0, ^°Ne, 
24Mg and 28Si. 

We finish by recalling the major problems that, among 
all the present intermediate-mass AGB models, remain to 
be solved. 

— We have not yet identified the ^^C additional source 
needed to quantitatively explain the s— process. This 
process is however required in order to account for 
many observations unquestionably indicating that 
AGB stars are responsible for the main component of 
the solar system heavy element distribution. Some sug- 
gestions have been made. Radiative diffusion operat- 
ing during the TP-AGB phase of low-mass stars could 
substantially help. Although already invoked, such a 
process has never been yet included in AGB star evo- 
lution models. 

— The rather high level of ^Li production detected at the 
surface of some evolved but relatively faint AGB stars 
of our galaxy still remains unaccountable. Thus, the 
Cameron-Fowler mechanism does not explain all the 
super-lithium-rich AGB stars. As these galactic ob- 
jects have relatively low initial masses, the explana- 
tion could also come from the radiative slow-transport 
processes above mentioned. 

— Last but not least, some C stars are observed, es- 
pecially in the galactic bulge, with low luminosities 
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Table 8. Total mass ejected (in Mq) in the wind for some nuclides at the RGB top, at the end of the E-AGB phase and at 
the AGB tip, just before the ejection of the planetary nebula. Also given is the net contribution for these elements, i.e. the 
total final mass ejected reduced by the amount of matter that would have been ejected with the initial composition. Here are 
presented elements from He to N 
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-7 


7.73 10"* 
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■5 
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■5 
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6.52 10"^ 




at AGB tip 
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3.03 10"" 
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-4 
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-5 


3.3010"" 


3.0410"" 


9.19 10"^ 
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2.25 10""* 


2.07 10"" 


3.46 10" 


.5 


-9.44 10" 




2.7410"" 
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-4.27 10"® 


^He 


at RG top 


6.5310"* 


4.28 10"" 


8.69 10" 


-4 


5.84 10" 


-4 


1.1510"" 


1.89 10"^ 
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6.91 10"^ 


9.59 10" 


-2 
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-2 


5.11 10"^ 


0.13 


0.18 




at AGB tip 


0.59 


0.91 


1.32 




1.59 




0.55 


0.72 


1.28 




net yield 


3.7610"^ 


3.20 10"^ 


0.17 




0.23 




3.5610"^ 


3.70 10"^ 


0.37 


^Li 


at RG top 


1.5910"" 


1.02 10"" 


2.49 10" 


11 


1.45 10" 


11 


1.4210"" 


2.6710"^^ 


1.2610"" 
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2.63 10"" 


3.3510"" 


5.18 10" 


11 


3.8710" 


11 


4.4710"^^ 


1.15 10"" 


1.96 10"" 




at AGB tip 


2.92 10"^° 


6.3510"" 


7.33 10" 


-8 


1.08 10" 


-7 


8.23 10"" 


4.6010"" 


3.5710"" 




net yield 
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-2.95 10"® 
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-8 
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-8 
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4.25 10"® 


9.11 10" 
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6 
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-4 
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-4 
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-3 
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-3 
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net yield 
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_2 
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-9 
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-9 
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2.90 10" 


-9 
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-9 
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"N 
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1.04 10"^ 
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1.9710" 


-6 
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at which the models fail at predicting the occurrence 
of the third dredge-up. This is probably due to our 
bad knowledge (and treatment) of finely tuned con- 
vection motions in boundaries of stellar regions having 
strong chemical composition gradients. Clearly, some 
kind of extra-mixing (deeper inside nuclearly processed 
regions) has to be invoked. However, no self-consistent 
mechanism has been settled and tested through com- 
plete evolution models up to now. 



Other confrontations with observations require the 
modeling of low-mass AGB stars. A clear distinction be- 
tween low- and intermediate-mass stars is justified. In- 
deed, numerous observations indicate that "non-standard" 
particle transport processes are acting inside low-mass 
stars at different phases of their evolution, that substan- 
tially modify the chemical structure compared to what is 
obtained in classical models. During the various dredge- 
up episodes, matter up-heaved to the surface has con- 
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Table 9. Same as Table 8 for elements from O to Na 
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3.96 10" 


-7 


3.8410" 


-7 


2«Ne 


at RG top 
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11 


8.38 10" 


12 
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-7 
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-8 
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-8 
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-8 
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-7 
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-7 
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-7 
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-7 
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-6 
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-7 
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-8 
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-7 
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-4 
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-5 
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-4 
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-4 
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-3 
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-4 
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-3 
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-5 


-1.51 10 
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1.11 10" 
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-4 
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7.76 10" 


-6 
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-4 
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-5 
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-5 
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-6 
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-6 
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-6 
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scqucntly a diffcaxmt composition. This is not the case 
in intermediate-mass stars. In conclusion, low-mass AGB 
stellar models, maintaining to make detailed predictions 
concerning the cvohition of surface isotopic ratios, have 
to include such slow-particle transport processes. They 
are currently being calculated and will be presented in 
a next-coming paper. 
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